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Abstract 

In this paper, we show the optimahty of a certain class of disturbance-afSne control policies in the 
context of one-dimensional, constrained, multi-stage robust optimization. Our results cover the finite 
horizon case, with minimax (worst-case) objective, and convex state costs plus linear control costs. We 
develop a new proof methodology, which explores the relationship between the geometrical properties 
of the feasible set of solutions and the structure of the objective function. Apart from providing an 
elegant and conceptually simple proof technique, the approach also entails very fast algorithms for the 
case of piecewise affine state costs, which we explore in connection with a classical inventory management 
application. 

1 Introduction. 

Multi-stage optimization problems under uncertainty have been prevalent in numerous fields of science and 
engineering, and have elicited interest from diverse research communities, on both a theoretical and a prac- 
tical level. Several solution approaches have been proposed, with various degrees of generality, tractability, 
and performance guarantees. Some of the most successful ones include exact and approximate dynamic 
programming, stochastic programming, sampling-based methods, and, more recently, robust and adaptive 
optimization, which is the focus of the present paper. 

The topics of robust optimization and robust control have been studied, under different names, by a 
variety of academic groups, mostly in operations research (Ben-Tal and Nemirovski [1999], Ben-Tal et al. 
[2002], Ben-Tal and Nemirovski [2002], Bertsimas and Sim [2003], Bertsimas and Sim [2004], Bertsimas 
et al. [2004]) and control theory (Bertsekas and Rhodes [1971], Fan et al. [1991], El-Ghaoui et al. [1998], 
Grieder et al. [2003], Bemporad et al. [2003], Kerrigan and Maciejowski [2004], Zhou and Doyle [1998], 
DuUerud and Paganini [2005]), with considerable effort put into justifying the assumptions and general 
modeling philosophy. As such, the goal of the current paper will not be to justify the use of robust (and, 
more generally, distribution-free) techniques. Rather, we will take the modeling approach as a given, and 
investigate tractability and performance issues in the context of a certain class of optimization problems. 
More precisely, we will be concerned with the following multi-stage decision problem: 

Problem 1.1. Consider the following one-dimensional, discrete, linear, time-varying dynamical system: 

Xk+i = ak ■ Xk -\- Pk ■ Uk -\- 7fc • Wk (1) 

where ak,Pk,lk are known scalars, and the initial state xi £ R is specified. The random disturbances 
Wk are unknown, but bounded: 

Wk eWk [Wk,Wk]. (2) 
We would like to find a sequence of robust controllers {uk}, obeying certain constraints: 

Uke[Lk,Uk] (3) 
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(Lk,Uk are known and fixed), and optimizing the following min-max cost function over a finite horizon 
1 T- 



lie/ 

;M = mm 

ui 



Ci • lii + max 



hi{x2) + • • • + min 



Ck ■ Uk + max [hk{xk+i) 



+min [ ct ■ Ut + max hT{xT+i) 
where the functions /i^ : R ^ R are known, convex and coercive, and Cfe > are fixed and known. 



(4) 



The problem corresponds to a situation in which, at every time step k, the decision maker has to compute 
a control action Uk, in such a way that certain constraints (3) are obeyed, and a cost penalizing both the state 
{hk{xk+i)) and the control {ck ■ Uk) is minimized. The uncertainty, Wk, always acts so as to maximize the 
costs, hence the problem solved by the decision maker corresponds to a worst-case scenario (a minimization of 
the maximum possible cost). Examples of such problems include the case of quadratic state costs, hk{xk+i) = 
rk+i ' {rk+i > 0), as well as norm-1 or norm-cx) costs, hk(xk+i) = rk+i ■ \xk+i \ , ('''fc+i > 0), all of which 
have been studied extensively in the literature in the unconstrained case (see, for example, Zhou and Doyle 
[1998], and DuUerud and Paganini [2005]). 

The solution to Problem 1.1 could be obtained using a "classical" Dynamic Programming formulation 
(see Bertsekas [2001]), in which the optimal policies ul{xk) and the optimal value functions J^ixk) are 
computed backwards in time, starting at the end of the planning horizon, k = T. The resulting policies are 
piecewise afSne in the state xt, and have properties that are well known and documented in the literature, 
dating back to Scarf et al. [1958]. 

In the current paper, we would like to study the performance of a new class of policies, where instead of 
regarding the controllers Uk as functions of the state Xk, one seeks direct paramcterizations in the observed 
disturbances: 

Wfc : Wi X W2 X • • • X Wk-i R. (5) 
In this framework, we would require that constraint (3) should be robustly feasible: 

Ukiw'')e[Lk,Uk], yw'^'^ {wi,...,Wk-i)eWiX---xWk-i. (6) 

Note that if we insisted on this category of paramcterizations, then we would have to consider a new state 
for the system, Xk, which would include all the past-observed disturbances, w''. Furthermore, while Xk 
summarizes all the past information needed to make an optimal decision at stage k, the same would not 
necessarily be true for w'^, so we may want to include even more information in Xk (for example, the 
previous controls {ut}i<t<fe or the previous states {a:;t}i<t<fe or some combination thereof). Compared with 
the original, compact state formulation, Xk, the new state Xk would become much larger, and solving the 
Dynamic Program with state variable Xk would produce exactly the same optimal objective function value 
JmM- Therefore, one would be tempted to ask what the benefit for introducing such a complicated state 
might be. 

The hope is that, by considering policies over a larger state, one might be able to obtain simpler functional 
forms, for instance, affine policies. These have a very compact representation, since only the coefficients of the 
parameterization are needed, and, for certain classes of convex costs hk{-), there may be efficient procedures 
available for computing them. 

This approach is also not new in the literature. It has been originally advocated in the context of stochastic 
programming (Rockafellar and Wets [1978]), then in robust optimization (Ben-Tal et al. [2004]), and extended 
to linear systems theory (Ben-Tal et al. [2005a], Ben-Tal et al. [2006]), with notable contributions from 
researchers in robust model predictive control and receding horizon control (see Bemporad et al. [2003], 
Kerrigan and Maciejowski [2004], Lofberg [2003], Skaf and Boyd [2008], and references therein). In all the 
papers, which usually deal with the more general case of multi-dimensional linear systems, the authors show 
how the reformulation can be done, and how the corresponding affine policies can be found by solving specific 
types of optimization problems, which vary from linear and quadratic programs (Ben-Tal et al. [2005a], 
Kerrigan and Maciejowski [2004], Kerrigan and Maciejowski [2003]) to conic and semi-definite (Ben-Tal 
et al. [2005a], Lofberg [2003], Bertsimas and Brown [2007], Grieder et al. [2003]), or even multi-parametric, 
linear or quadratic programs (Bemporad et al. [2003]). The first steps towards analyzing the properties of 
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such parameterizations were made in Kerrigan and Maciejowski [2004], where the authors show that, under 
suitable conditions, the resulting affine parameterization has certain desirable system theoretic properties 
(stability and robust invariancc). Another notable contribution was made by Goulart and Kerrigan [2005], 
who prove that the class of affine disturbance feedback policies is equivalent to the class of affine state 
feedback policies with memory of prior states, thus subsuming the well known classes of open- loop and pre- 
stabilizing control policies. However, to the best of our knowledge, apart from these theoretical advances, 
there has been very little progress in proving results about the quality of the objective function values 
resulting from the use of such parameterizations. 

Our main result, summarized in Theorem 3.1 of Section 3, is that, for Problem 1.1 stated above, affine 
policies of the form (5) are, in fact, optimal. Furthermore, we are able to prove that a certain (affine) 
relaxation of the state costs is also possible, without any loss of optimality, which gives rise to very efficient 
algorithms for computing the optimal affine policies when the state costs /ifc(-) are piece- wise affine. To the 
best of our knowledge, this is the first result of its kind, and it provides intuition and motivation for the 
widespread advocation of such policies in both theory and applications. Our theoretical results are tight 
(if the conditions in Problem 1.1 are slightly perturbed, then simple counterexamples for Theorem 3.1 can 
be found), and the proof of the theorem itself is atypical, consisting of a forward induction and making 
use of polyhedral geometry to construct the optimal affine policies. Thus, we are able to gain insight into 
the structure and properties of these policies, which we explore in connection with a classical inventory 
management problem. 

The paper is organized as follows. Section 2 presents an overview of the Dynamic Programming formula- 
tion in state variable Xk, extracting the optimal policies M^(xfc) and optimal value functions J^(xk), as well as 
some of their properties. Section 3 contains our main result, and briefly discusses some immediate extensions 
and computational implications. In Section 4, we introduce the constructive proof for building the affine 
control policies and the affine cost relaxations, and present counterexamples that prevent a generalization 
of the results. In Section 5, we explore our results in connection with a classical inventory management 
problem. Section 6 presents our conclusions and directions for future research. 



2 Dynamic Programming Solution. 

As already mentioned in the introduction, the solution to Problem 1.1 can be obtained using a "classical" 
Dynamic Programming formulation [Bcrtsekas, 2001], in which the state is taken to be Xk-, and the optimal 
policies u*f,{xk) and optimal value functions J^{xk) are computed starting at the end of the planning horizon, 
k = T, and moving backwards in time. In this section, we will briefly outline the solution of the Dynamic 
Program for our problem, and will state some of the key properties that will be used throughout the rest of 
the paper. For completeness, a full proof of the results is included in Section 7.1 of the Appendix. 

In order to simplify the notation, we remark that, since the constraints on the controls Uk and the 
bounds on the disturbances Wk are time- varying, and independent for different time-periods, we can restrict 
attention, without loss of generality^, to a system with ak = Pk — Jk — With this simpliflcation, the 
problem that we would like to solve is the following: 



(DP) 



JmM = min 

"1 



ci ■ ui + max 



hiix2) ^ h min 



s.t. Xk+i = Xk + Uk + Wk 
Lk <Uk < Uk 
Wk G Wk = [Wfc ,Wk]. 



Ck-Uk + max [hk{xk+i) + ■ . 



-min [ ct ■ Ut + max hT{xT+i] 

Ut \ Wt 



VfcG{i,2,...,r} 



^Such a system can always be obtained by the linear change of variables ij. = — — , and by suitably scaling the bounds 

_ ni=i 
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The corresponding Bellman recursion for (DP) can then be written as follows: 



mm 

Lk<uk<Uk 



where J^_^_i{xt+i) = 0. By defining: 



Ck-Uk+ max [hkixk+Uk + Wk) + Jk+i{xk+Uh + Wk) 



Vk = Xk + Uk 
gk iUk) = max [ hkiuk + Wk) + Jk+i iVk + Wk) 



(7) 



we obtain the following solution to the Bellman recursion (see Section 7.1 in the Appendix for the derivation): 



ul{xk) 



Uk, if Xk <y_^ - Uk 

~Xk + y*, otherwise 
Lk, if Xk >yk-Lk 

Ck ■ Uk + gk{xk + Uk), if Xk<y^ - Uk 



Jlixk) = Ck ■ ul{xk) + gk {xk + ul{xk)) = {ck- {y* -Xk)+ gk{y*), otherwise 

Ck-Lk+gkixk+Lk), if Xk>yk~ Lk 



(8) 



(9) 



where y* G [y^ , y^,], and [y^ , ?/^,] represents the (compact) set of minimizers of the convex function Ck-y+gkiy)- 
A typical example of the optimal control law and the optimal value function is shown in Figure 1. 




y*-Uk y*-Lk 



Xk 




y*-Uk y*-L| 



Figure 1: Optimal control law ul{xk) and optimal value function J^{xk) at time k. 



The main properties of the solution that will be relevant for our later treatment are listed below: 
PI The optimal control law u'^(xk) is piecewise afHne, with 3 pieces, continuous and non-increasing. 
P2 The optimal value function, J^(xk), and the function gk{yk) are convex. 

P3 The difference in the values of the optimal control law at two distinct arguments s <t always satisfies: 
- "fc(0 = -fk -{s-t), for some fk G [0, 1]. 

P4 The optimal value function, J^{xk), has a subgradient at most —Ck for Xk < y^, — Uk-, exactly —Ck in 
the interval [y^ — Uk,yk ~ Lk], and at least —Ck for Xk > Ijk — Lk- 

P5 The function gk{yk) is decreasing, with a subgradient at most —Ck, in the interval {—oo,y^, and is 
increasing, with a sugradient at least —Ck, in the interval [y^ , oo). 
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3 Optimality of AfRne Policies in w^. 

In this section, we introduce our main contribution, namely a proof that policies that are affine in the 
disturbances w'' are, in fact, optimal for problem {DP). Using the same notation as in Section 2, we can 
summarize our main result in the following theorem: 

Theorem 3.1. For every time step k ~ 1, . . . ,T, the following quantities exist: 
• An ajfine control policy: 

k-l 

qk{wi^) = qkfi + ^ qk,t ■ wt (10) 



An affine running cost: 



t=i 



(11) 



such that the following properties are obeyed: 
Lk < Qkiw'') < Uk, 



Zk{w''+^) > hk ixi+Y,{9t{w') 



Wt 



Vw'' G Wi X •■• X Wk-i 
Vio'^+i e Wi X ••■ X Wfc 



J, 



niM 



max 

U)l,...,lUi. 



Wt 



(12) 
(13) 

(14) 



Let us interpret the main statements and results in the theorem. Equation (12) confirms the existence 
of an affine policy qk{w^) that is robustly feasible, i.e. that obeys the control constraints, no matter what 
the realization of the disturbances may be. Equation (13) states the existence of an affine cost Zk{w^'^^) 
that is always larger than the convex state cost hk{xk+i) incurred when the affine policies {qt{')}i<t<k are 
used. Equation (14) guarantees that, despite using the (suboptimal) affine control law qk{w^)^ and incurring 
a (potentially larger) affine stage cost Zk{w^'^^), the overall objective function value JmM is, in fact, not 
increased. This translates in the following two main results: 

• Existential result. Affine policies qk{w^) are, in fact, optimal for Problem 1.1. 

• Computational result. When the convex costs hk{xk+i) are piecewise affine, the optimal affine policies 
{^qk{w^)^ ^^i^^rp can be computed by solving a Linear Programming problem. 

To see why the second implication would hold, suppose that hk{xk+i) is the maximum of nik affine functions: 

hk{xk+i) = max (jp\. ■ Xk+i +pLo) ' « G {1, ■ • -^mk}. 

Then an optimal affine policy qk{w^) can be obtained by solving the following optimization problem (see 
Ben-Tal et al. [2005a,b]): 

min J 

J;{qkd}]{zk,t} 

s.t. Vu; G Wi X W2 X ••• X Wt, V/c = l,...,r: 



{AARC) 

Vi G {1,, 



fc=i 



fc-i 



Cfc • qkX) + Zkfi + '^{cf qk,t + Zk,t) ■ Wt + Zk,k ■ Wk 



,mk} 



Zk 



,0 + E ^k,t ■ Wt > pI 



Xi 



E 



9t,o 



E 9*'^ ' ^■^ 



■ Wt 



- Pk,o 



fe-i 



Lk < Qkfi + E ■Wt<Uk. 
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Although Problem {AARC) is still a semi-infinite LP (due to the requirement of robust constraint feasibility, 
Vio), since all the constraints are inequalities that are bi-affine in the decision variables and the uncertain 
quantities, a very compact reformulation of the problem is available. In particular, with a typical constraint 
in (AARC) written as: 

T 

Ao[a;] + ^ At[a;] • wt < 0, Vi(? £ Wi x W2 x • • • x Wt , 

t=i 

where X[x] denotes that A is an affine function of the decision variables x, it can be shown (see Ben-Tal and 
Nemirovski [2002], Ben-Tal et al. [2004] for details) that the previous condition is equivalent to: 



'Ao[^] + ELi (a. W . ^ -1- ^ . 6) < 
~^t<Mx]<^t, t=l,...,T, 



(15) 



which are linear constraints in the decision variables x,^. Therefore, {AARC) can be reformulated as a 
Linear Program, with O {T"^ ■ maxfc mk) variables and O (T^ • max^ mk) constraints, which can be solved 
very efficiently using commercially available software. 

We will conclude our observations by making one last remark related to an immediate extension of the 
results. Note that in the statement of Problem 1.1, there was no mention about constraints on the states Xk 
of the dynamical system. In particular, one may want to incorporate lower or upper bounds on the states, 
as well: 

Ll<xk< (16) 

We claim that, in case the mathematical problem including such constraints remains feasible^, then affine 
policies will, again, be optimal. The reason is that such constraints can always be simulated in our current 
framework, by adding suitable convex barriers to the stage costs hk{xk+i)- In particular, by considering the 
modified, convex'^ stage costs: 

hk{xk+i) = hk{xk+i) + l[Li^^,u-^^]{xk+i), where 

1 c/x (xk+i) = 1°' ""'^^^ ^ '^^+^' ^^"+^^ 
[-f-fc+i.t^fe+il '^"+1 [00, otherwise, 

it can be easily seen that the original problem, with convex stage costs hk{-) and state constraints (16), 
is equivalent to a problem with the modified stage costs hk{-) and no state constraints. And, since affine 
policies are optimal for the latter problem, the result is immediate. Therefore, our decision to exclude such 
constraints from the original formulation was made only for sake of brevity and conciseness of the proofs, 
but without loss of generality. 



4 Proof of Main Theorem. 

The current section contains the proof of Theorem 3.1. Before presenting the details, we would first like to 
give some intuition behind the strategy of the proof, and introduce the organization of the material. 

Unlike most Dynamic Programming proofs, which utilize backward induction on the time-periods, we will 
have to proceed with a forward induction. Section 4.1 will present a test of the first step of the induction, 
and will then introduce a detailed analysis of the consequences of the induction hypothesis. 

We will then separate the completion of the induction step into two parts. In the first part, handled in 
Section 4.2, by exploiting the structure provided by the forward induction hypothesis, and making critical 
use of the properties of the optimal control law u*^.{xk) and optimal value function J^ixk) (the solutions 
to problem {DP) in Section 2), we will be able to introduce a candidate affine policy qk{w^). In Section 

^Such constraints may lead to infeasible problems. For example: T = l,x\ = 0,u\ € [0,1], toi G [0,1], X2 £ [5,10]. 
^Thc functions /ife(-) arc convex in since l{[Lfe^i, l/j-^i]), the indicator function of a convex set, is convex (see pages 

28, 33 in Rockafellar [1970] or Example 7.1.2 in Bertsekas [2003]), and the sum of convex functions remains convex. 
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4.2.1, we will then prove that this policy is robustly feasible, and preserves the min-max value of the overall 
problem, JmM, when used in conjunction with the original, convex state costs, hk{xk+i)- 

Similarly, for the second part of the inductive step, (Section 4.3), by re-analyzing the feasible sets of the 
optimization problems resulting after the use of the (newly computed) afBne policy qk{w''), we will determine 
a candidate affine cost Zk{w^~^^), which wc will prove (Section 4.3.1) to be always larger than the original 
convex state costs, hk{xk+i)- However, despite this fact, in Section 4.3.1 wc will also show that when this 
affine cost is incurred, the overall min-max value JmM remains unchanged, which will complete the proof of 
the inductive step. 

Section 4.3.2 will then conclude the proof of Theorem 3.1, and will outline several counterexamples that 
prevent an immediate extension of the result to more general cases. 



4.1 Induction Hypothesis. 

As mentioned before, the proof of the theorem will utilize a forward induction on the time-step fc. Let us 
begin by verifying the induction at fc = 1. 

Using the same notation as in Section 2, by taking the afRne control to be qi w|(xi), we immediately 
get that qi, which is simply a constant, is robustly feasible, so (12) is obeyed. 

Furthermore, since qi = ul{xi) is optimal, we can write the overall optimal objective value, JmM, as: 

J7nM = Ji{^i)= ™n [ci • Ml -I- gi(a::i ui) ] = ci • gi +51 (xi + gi) 

uie[Li,Ui] 

ci • gi -I- max{ (ft-i + J2 ) (a;i -I- gi + Wi) , {hi + J2) {xi + qi +wi)}. (17) 

Next, we introduce the affine cost zi(u;i) == zi.o + zi^i ■ wi, where we constrain the coefficients zi^i to satisfy 
the following system of linear equations: 

\ zifi + zi^i ■ = hi{xi+qi+w^) 
[zifi + zi^i -Wi = hi{xi + qi +Wi) . 

Note that for fixed xi and qi, the function zi{wi) is nothing but a linear interpolation of the mapping 
wi I— > hi{xi + qi + wi), matching the value at points {wi ,wi}. Since hi is convex, the linear interpolation 
defined above clearly dominates it, so condition (13) is readily satisfied. Furthermore, by (17), JmM is 
achieved for wi S {w^ ^wi}, so condition (14) is also obeyed. 

Having checked the induction at time fc = 1, let us now assume that the statements of Theorem 3.1 are 
true for times t — 1, . . . , fc. Equation (14) written for stage fc then yields: 



'mM = max 

Wl,...,Wk 



k / k \ 

(ct ■ qt{w*) + zt {w'+^)) + Ui + ^ {qt{w') + Wt) 

t=i \ t=i J 

= max [0i-fJ,!+i((?2)], (18) 

{k k \ 

(01, 02) : Bi ^ (ct • qtiw') {w'+^)) , 02 = ^1 + E Wt)\ . (19) 

t=\ t=l J 

Since {'Zt}i<t<fe and {•Zt}i<t<j; arc affine functions, this implies that, although the uncertainties lui, . . . , Wfe 
lie in a set with 2*^ vertices (the hyper-rectangle Wi x • • • x Wfc), they are only able to affect the objective JmM 
through affine combinations taking values in the set 0. Such a polyhedron, arising as a 2-dimensional affine 
projection of a fc-dimensional hyper-rectangle, is called a zonogon. It belongs to a larger class of polytopes, 
known as zonotopes, whose combinatorial structure and properties are well documented in the discrete and 
computational geometry literature. The interested reader is referred to Chapter 7 of Ziegler [2003] for a very 
nice and accessible introduction. 

The main properties of a zonogon that we are interested in are summarized in Lemma 7.2, found in the 
Appendix. In particular, the set Q is centrally symmetric, and has at most 2fc vertices (see Figure 2 for an 
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Figure 2: Zonogon obtained from projecting a hypercube in R". 



example). Furthermore, by numbering the vertices of in counter-clockwise fashion, starting at: 

Vq = = argmax< argmin{0 G 0} > , (20) 
ei I 02 J 

we can establish the following result concerning the points of that are relevant in our problem: 

Lemma 4.1. The maximum value in optimization problem (18) is achieved for (61,92) {vo,vi, ... ,Vk}. 

Proof. The optimization problem described in (18) and (19) is a maximization of a convex function over 
a convex set. Therefore (see Section 32 of Rockafellar [1970]), the maximum is achieved at the extreme 
points of the set O, namely on the set {vo,Vi, . . . ,V2p-i,V2p = vq}, where 2p is the number of vertices of 
0. Letting O denote the center of 0, by part (iii) of Lemma 7.2 in the Appendix, we have that the vertex 
symmetrically opposed to v~ , namely = 20 — , satisfies = Vp. 

Consider any vertex Vj with j G {p + I, . . . ,2p ~ 1}. From the definition of , , for any such vertex, 
there exists a point 7;^ G [v^ , v^], with the same 6'2-coordinate as Vj, but with a (?i-coordinate larger than 
Vj (refer to Figure 2). Since such a point will have an objective in problem (18) at least as large as vj, 
and v'j' G [i;o,Wp], we can immediately conclude that the maximum of problem (18) is achieved on the set 
{vq, . . . , Vp}. Since 2p < 2k (see part (ii) of Lemma 7.2), we can immediately arrive at the conclusion of the 
lemma. □ 

Since the argument presented in the lemma will be recurring throughout several of our proofs and 
constructions, we will end this subsection by introducing two useful definitions, and generalizing the previous 
result. 

Consider the system of coordinates (6*1,6*2) in R^, and let S C denote an arbitrary, finite set of 
points and V denote any (possibly non-convex) polygon such that its set of vertices is exactly S. With 
y~ = argmaxg^ jargmiug^ {0 G 5}} , = argmaxg^ jargmaxg^ {0 G 5}}, by numbering the vertices of 

the convex hull of S in a counter-clockwise fashion, starting at j/q = y^, and with — y+, we define the 
right side of V and the zonogon hull of S as follows: 

Definition 4.1. The right side of an arbitrary polygon V is: 

r-side(7') ='{yo,2/i,.-.,y™}. (21) 
Definition 4.2. The zonogon hull of a set of points S is: 

z-huU (5) =' |y G R2 : y = y^ + Y,w,- (y,^ - y^_^) , < < 1 1 . (22) 
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-Ti ' ' BT 

Figure 3: Examples of zonogon hulls for different sets S G R^. 



Intuitively, r-side(T') represents exactly what the names hints at, i.e. the vertices found on the right side 
of v. An equivalent definition using more familiar operators could be the following: 



r-side(7') = ext cone 





-1 



■{V) 



where cone(-) and conv(-) represent the conic and convex hull, respectively, and ext(-) denotes the set of 
extreme points. 

Using Definition 7.1 in Section 7.2 of the Appendix, one can see that the zonogon hull of a set S 
is simply a zonogon that has exactly the same vertices on the right side as the convex hull of S, i.e. 
r-side (z-huU (iS)) = r-side (conv (5)). Some examples of zonogon hulls are shown in Figure 3 (note that the 
initial points in S do not necessarily fall inside the zonogon hull, and, as such, there is no general inclusion 
relation between the zonogon hull and the convex hull). The reason for introducing this object is that it 
allows for the following immediate generalization of Lemma 4.1: 

Corollary 4.1. IfV is any polygon in = 62) with a finite set S of vertices, and / : R ^ R is any 
convex function, then the following string of equalities holds: 



max I oi 



Oi + f{02)\ 



max \0i + f(e2)]= max \ei + f(02)]= max [^i 
(6(i,6(2)econv(p) (ei.e2)es ' (ei,e2)er-sidcCP) 



fi02)\ 



U{92)]= max [0i+/(02)" 

(e'i,e'2)er-sido(z-hull(5)) 

Proof. The proof is identical to that of Lemma 4.1, and is omitted for brevity. 



max 

(6li,6(2)ez-hull(5) 



□ 



Using this result, whenever we will be faced with a maximization of a convex function 0i + /(02), we will be 
able to switch between different feasible sets, without affecting the overall optimal value of the optimization 
problem. 

In the context of Lemma 4.1, the above result allows us to restrict attention from a potentially large set 
of relevant points (the 2'^ vertices of the hyper- rectangle Wi x • • • x Wk), to the fc + 1 vertices found on 
the right side of the zonogon Q, which also gives insight into why the construction of an affine controller 
qk+ii'w'''^^) with k + 1 degrees of freedom, yielding the same overall objective function value JmM, might 
actually be possible. 

In the remaining part of Section 4.1, we would like to further narrow down this set of relevant points, 
by using the structure and properties of the optimal control law w^^.^(xfc+i) and optimal value function 
Jk+ii^'k+i)} derived in Section 2. Before proceeding, however, we will first reduce the notational clutter by 
introducing several simplifications and assumptions. 

4.1.1 Simplified Notation and Assumptions. 

To start, we will omit writing the time subscripts (k or k + 1) whenever possible, and will take: 

6'i(ti;) = ao + • w,;; 6*2(10) = &o + ^ • Wi; qk+i{w) = q{w) = qq +^^qi ■ Wi , (23) 



i=l 



4=1 
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where a ~ (oi, . . . , a^.) and b = (61, ... , bk) are the generators of the zonogon O. We will use the same 
counter-clockwise numbering of the vertices of O as introduced earlier in Section 4.1: 

t>0 = V^,...,Vp = V+, . . . ,t!2p = v~ , (24) 

where 2p is the number of vertices of 0, and we will use the following overloaded notation for 61 and 92'. 

• 9i{w), for w E R'^, will denote the 9i value assigned by the affinc projection to w, via (23). 

• 9i [v], applied to v G R^, will be used to denote the 6'i-coordinate of the point v. 

Also, since 02 = Xk+i, instead of referring to J^^^{xk+i) and ul^^{xk+i), we will use J*(02) and u*{92), and 
we will occasionally use the short-hand notations u*{vi). J*{vi) and g (vi), instead of u* {92[vi]), J* {92[vi]) 
and g {92[vi] + u*{92[vi])), respectively. 

Since many of the sets of interest will lie in R^, for a system with coordinates {x,y), we define the 
following convenient notation for the cotangent of the angle formed by an oriented line segment [M, N] with 
the X-axis: 

cotan (M, N) = ^EIL^I^^ where M = (xm, 2/j\/) e R^ N= (xjv, yjv) e R^- (25) 

We will also make the following simplifying assumptions: 

Assumption 1. The uncertainty vector at time k, = (wi, . . . ,Wk), belongs to the unit hypercube Tiu of 
R'', i.e. Wi = 0,Wi = 1, V i = 1, . . . , k. 

Assumption 2. The zonogon has a maximal number of vertices, i.e. p = k. 

Assumption 3. The vertex of the hypercube projecting to Vi, i £ {0, . . . , k}, is exactly [1, 1, . . . , 1, 0, . . . , 0], 
i.e. 1 in the first k components and thereafter (see Figure 2). 

These assumptions are made only to facilitate the exposition, and result in no loss of generality. To 
see this, note that the conditions of Assumption 1 can always be achieved by adequate translation and 
scaling of the generators a and b (refer to Section 7.2 of the Appendix for more details), and Assumption 
3 can be satisfied by renumbering the coordinates of the hyper-rectangle, i.e. renumbering the disturbances 
wi, . . . ,Wk. As for Assumption 2, we argue that an extension of our construction to the degenerate case 
p < fc is immediate (one could also remove the degeneracy by applying an infinitesimal perturbation to the 
generators a or b, with infinitesimal cost implications). 



4.1.2 Further Analysis of the Induction Hypothesis. 

With this simplified notation, by using (9) to express J*(-) as a function of u*{-) and g{-), we rewrite equation 
(18) as follows: 



(OPT) JrnM = max [71+5(72) ] , 
(7i.72)ef 

where f = | (71, 72) : 71 = ^^i + c • m*(02), 72 = ^2 + u*(02), (0i,02)ee} 



A characterization for the set F = (71,72) can be obtained by replacing the optimal, piecewise affine control 
law u*{92), given by^ (8) in equation (26): 

{9i+c- U, 92 + U) , if92<y* -U 

(71,72) ="( (6*1 - c • 6*2 + c ■ y*, y*) , otherwise (27) 
{9i+c- L, 92 + L) , ii92>y*~L 

Using the same overloaded notation^ for 71,72: 

*For simplicity, wc focus on the case when g{-) has a unique minimizer, such that y =y = y* in (8), (9). The analysis could 
be immediately extended to the set of minimizers, [3/ 

^To avoid confusion, in addition to different types of parentheses, we will consistenly use tilded arguments to denote points 
Vi G r, and non-tilded arguments for -u G 0. 
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• 'Yiiv), applied to t> G 8 C R^, will refer to the 7^ value assigned by the mapping (27) to v 

• 7,; [v], applied to i; G r C R^, will denote the 7,; coordinate of the point i), 

we can now provide a compact characterization for the maximizers in problem (OPT) from (26): 
Lemma 4.2. The maximum in problem (OPT) over T is reached on the right side of: 

Ar ='conv({i}o,...,i>fc}) , (28) 

where: 

ir = i^iiv,), Uvr) ) = idi[v,]+c-u*iv,), 92[v,] + u* (v,) ) , Z G {0, . . . , fc}. (29) 

Proof. By Lemma 4.1, the maximum in (18) is reached at one of the vertices vo,vi, . . . ,Vk of the zonogon 
O. Since this problem is equivalent to problem (OPT) in (26), written over F, we can immediately conclude 
that the maximum of the latter problem is reached at the points {'5i}i<i<fc given by (29). Furthermore, 
since g(-) is convex (see Property P2 of the optimal DP solution, in Section 2), we can apply Corollary 4.1, 
and replace the points Vi with the right side of their convex hull, r-side (Ar), without changing the result of 
the optimization problem, which completes the proof. □ 

Since this result will be central to our future construction and proof, we will spend the remaining part 
of the subsection discussing some of the properties of the main object of interest, the set of points on the 
right side of Ar, r-side(Ar). To understand the geometry of the set Ar and its connection with the optimal 
control law, note that the mapping (27) from to F will discriminate points 6 ^ (6*1, 6*2) £ 9 depending on 
their position relative to the horizontal band: 

Blu - { (^^1,^2) G R2 : 02 e [y* -U,y*-L]}. (30) 

In particular, from (27) and the definition of Vq, • . • , ffc in (24) and (20), we can distinguish a total of 
four distinct cases. The first three, shown in Figure 4, are very easy to analyze: 




Si 01 di 



Figure 4: Trivial cases, when zonogon lies entirely [CI] below, [C2] inside, or [C3] above the band Blu- 



[CI] If the entire zonogon Q falls below the band Blu, i-c 02 [vk] < V* — U, then f is simply a translation 
of 6, by (c • U, U), so that r-side (Ar) = {vq, Vi, . . . , Vk}- 

[C2] If 6 lies inside the band BlUi i-e. y* - U < 62 [vq] < 62 [vk] < y* — L, then aU the points in f wiU 
have 72 = y* , so F will be a line segment, and |r-side(Ar)| = 1. 

[C3] If the entire zonogon Q falls above the band Blu, i-e. O2 [vq] > y* — L, then F is again a translation 
of 6, by (c • L, L), so again r-side (Ar) = {vq, Vi, . . . , Vk}- 

The remaining case, [C4], is when O intersects the horizontal band Blu in nontrivial fashion. We can 
separate this situation in the three sub-cases shown in Figure 5, depending on the position of the vertex 
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71 71 71 



Figure 5: Case [C4]. Original zonogon B (first row) and tiie set F (second row) wlien Vt falls (a) under, (b) inside 
or (c) above the band Blu- 



Vt G r-sidc(0), where the index t relates the per-unit control cost, c, with the geometrical properties of the 
zonogon: 

. f f , if < c 

I max |ze{l,...,A:}:^>c| , otherwise . 

We remark that the definition of t is consistent, since, by the simplifying Assumption 3, the generators a, b 
of the zonogon Q always satisfy: 

> >■■■> (32) 
bi,b2,...,bk>0. 

An equivalent characterization of Vt can be obtained as the result of an optimization problem: 

Vt = argmin < argmax {9i — c ■ 62} > ■ (33) 
S2 [{0i,e2)ee J 

The following lemma summarizes all the relevant geometrical properties corresponding to this case: 

Lemma 4.3. When the zonogon Q has a non-trivial intersection with the band Bm (case [C4]), the convex 
polygon Ar and the set of points on its right side, r-side(Ar), satisfy the following properties: 

1. r-side(Ar) is the union of two sequences of consecutive vertices (one starting at Vq, and one ending at 
Vk), and possibly an additional vertex, Vt'- 

r-side(Ar) = {vo, Vi,...,Vs}iJ {vt} U {vr,Vr+i ...,vu}, for some s < r e {Q, . . . ,k] . 

2. With cotan(-, •) given by (25) applied to the (71,72) coordinates, we have that: 

J cotan [vs, i>min(t,r)) > ff^ , whenever t > s 



[cotan ({>max(t,s), «r) < f^, whenever t < 



(34) 



12 



While the proof of the lemma is slightly technical, which is why we have decided to leave it for Section 
7.3 of the Appendix, its implications are more straightforward. In conjuction with Lemma 4.2, it provides 
a compact characterization of the points Vi G T which are potential maximizcrs of problem (OPT) in (26), 
which immediately narrows the set of relevant points Vi G O in optimization problem (18), and, implicitly, 
the set of disturbances w G Tik that can achieve the min-max cost JmM- 

4.2 Construction of the AfRne Control Law. 

Having analyzed the consequences that result from using the induction hypothesis of Theorem 3.1, we now 
return to the task of completing the inductive proof, which amounts to constructing an affinc control law 
qk+i{w''^^) and an affine cost Zk+i{w''~^'^) that verify conditions (12), (13) and (14) in Theorem 3.1. We will 
separate this task into two parts. In the current section, we will exhibit an affine control law qk+i{'w'^^^) that 
is robustly feasible, i.e. satisfies constraint (12), and that leaves the overall min-max cost JmM unchanged, 
when used at time fc + 1 in conjunction with the original convex state cost, hk+iixk+2)- The second part of 
the induction, i.e. the construction of the affine costs Zfc+i(io'^+^) satisfying (13) and (14). will be left for 
Section 4.3. 

In the simplified notation introduced earlier, the problem we would like to solve is to find an affinc control 
law q{w) such that: 

JmM ^ max [ Oi + c ■ q{w) + g {02 + q{w)) ] 
we[as]'' 

L<q{w)<U, yweUk- 

The maximization represents the problem solved by the disturbances, when the affine controller, q{w), 
is used instead of the optimal controller, u*{92). As such, the first equation amounts to ensuring that the 
overall objective function JmM remains unchanged, and the inequalities are a restatement of the robust 
feasibility condition. The system can be immediately rewritten as: 

{AFF) JmM ^ max [71+5(72)] 

(7l,72)6r ^35^ 

L < q{w) < U , 

where: 

r = { (71, 72) : 71 - di+c- q{w), 72 - 62 + q{w), {61,62) £ o} . (36) 

With this reformulation, all our decision variables, i.e. the affine coefficients of q{w), have been moved 
to the feasible set F of the maximization problem {AFF) in (35). Note that, with an affine controller 
q{w) = qa + q^w, and 61,62 affine in w (23), the feasible set F will represent a new zonogon in R^, with 
generators given by a + c ■ q and b + q. Furthermore, since the function g is convex, the optimization 
problem {AFF) over F is of the exact same nature as that in (18), defined over the zonogon Q. Thus, in 
perfect analogy with our discussion in Section 4.1 (Lemma 4.1 and Corollary 4.1), we can conclude that the 
maximum in {AFF) will occur at a vertex of F found in r-side(F). 

In a different sense, note that optimization problem {AFF) is also very similar to problem {OPT) in 
(26), which was the problem solved by the uncertainties w when the optimal control law, u*{62), was used at 
time fc + 1. Since the optimal value of the latter problem is exactly equal to the overall min-max value, JmM, 
we will interpret the equation in (35) as comparing the optimal values in the two optimization problems, 
{AFF) and {OPT). 

As such, note that the same convex objective fimction, ^1 -I- 5(^2), is maximized in both problems, but 
over different feasible sets, F for {OPT) and F for {AFF), respectively. From Lemma 4.2 in Section 4.1.2, 
the maximum of problem {OPT) is reached on the set r-side(Ar), where Ar = coiiv ({i)o, Vl^ ■ ■ ■ , Vk})- From 
the discussion in the previous paragraph, the maximum in problem {AFF) occurs on r-sidc(F). Therefore, 
in order to compare the two results of the maximization problems, we must relate the sets r-side(Ar) and 
r-sidc(F). 

In this context, we introduce the central idea behind the construction of the affine control law, q{w). 
Recalling the concept of a zonogon hull introduced in Definition 4.2, we argue that, if the affine coefficients of 
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the controller, qo, q, were computed in such a way that the zonogon T actually corresponded to the zonogon 
hull of the set {vo,Vi, . . . ,Vk}, then, by using the result in Corollary 4.1, we could immediately conclude 
that the optimal values in (OPT) and {AFF) are the same. 

To this end, we introduce the following procedure for computing the afhnc control law q{w): 



Algorithm 1 Compute afhne controller q{w) 
Require: ei{w),e2{w), g{-),u*{-) 

1: if (6 falls below Blu) or (9 C Blu) or (9 falls above Blu) then 

2: Return q{w) ^ u*{92{w)). 

3; else 

4; Apply the mapping (27) to obtain the points Vi, z G {0, . . . , k}. 
5: Compute the set Ar = conv {{vo, ■ ■ ■ , Vk})- 

6: Let r-side(Ar) = {vq, vi, . . . , Vg} U {vt} U {vr, ■ ■ ■ , Vk} be the set of points on the right side of Ap. 
7: Solve the following system for qo, . . . ,qk and Ku, K^: 

Vt>i G r-side(Ar) (matching) 
Vz G {s + 1, . . . , min(t, r)} (alignment below t) ^g, 

Vi G {max(t, s) + 1, . . . , ?'} (alignment above t) 

8: Return q{w) = qn + J2i^i qiWi. 
9: end if 



' go H h <7i = u* {vi) 

a., + c- q, 



bi + qi 
a, + c ■ 

bi + a: 



= Ku, 



Before proving that the construction is well-defined and does produce the expected result, we will give 
some intuition for the constraints in system (37). In order to have the zonogon T be the same as the zonogon 
hull of {vo, ■ ■ ■ ,Vk}, must ensure that the vertices on the right side of F exactly correspond to the points on 
the right side of Ar = conv {{vo, . . . , Vk})- This is achieved in two stages. First, we ensure that vertices Wi 
of the hypercube Tik that are mapped by the optimal control law u*{-) into points Vi G r-side(Ar) (through 

(19) (29) 

the succession of mappings Wi ^ Vi £ r-side(9) ^ Vi £ r-side(Ar)), will be mapped by the affine control 

law, q{wi), into the same point Vi (through the mappings Wi Vi G r-side(9) Vi G r-side(Ar)). 
This is done in the first set of constraints, by matching the value of the optimal control law at any such 
points. Second, we ensure that any such matched points Vi actually correspond to the vertices on the right 
side of the zonogon F. This is done in the second and third set of constraints in (37), by computing the 
affine coefficients qj in such a way that the resulting segments in the generators of the zonogon F, namely 

'^^ I , are all aligned, i.e. have the same cotangent, given by the Ku, Kl variables. Geometrically, 
~ 1j ) 

this exactly corresponds to the situation shown in Figure 6 below. 

We remark that the above algorithm does not explicitly require that the control q{w) be robustly feasible, 
i.e. second condition in (35). However, we will soon prove that this is a direct result of the way matching 
and alignment are performed in Algorithm 1. 



4.2.1 Affine Controller Preserves Overall Objective and Is Robust. 

In this section, we prove that the affine control law q(w) produced by Algorithm 1 satisfies the requirements of 
(35), i.e. it is robustly feasible, and it preserves the overall objective function JmM when used in conjunction 
with the original convex state costs, /i(-). With the exception of Corollary 4.1, all the key results that we will 
be using are contained in Section 4.1.2 (Lemmas 4.2 and 4.3). Therefore, we will preserve the same notation 
and case discussion as initially introduced there. 

First consider the condition on line 1 of Algorithm 1, and note that this corresponds to the three trivial 
cases [CI], [C2] and [C3] of Section 4.1.2. In particular, since 6*2 = Xfe+i, we can use (8) to conclude that 
in these cases, the optimal control law «*(•) is actually affine: 
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Vo = V 



Original zonotopc 6. (9i Set f and points yj e r-sidc(r). 7i 

Figure 6: Outcomes from the matching and alignment performed in Algorithm 1. 

[CI] If falls below the band BlUi then the upper bound constraint on the control at time k is always 
active, i.e. the optimal control is u*{92{w)) ^ U, Vtu G Hk- 

[C2] If C Blu, then the constraints on the control at time k are never active, i.e. u*{92{w)) — y* — 6*2(11'), 
hence afRne in w, since 62 is affine in w, by (23). 

[C3] If falls above the band Blu, then the lower bound constraint on the control is always active, i.e. 
u*{92{w)) = L, yweHk- 

Therefore, with the assignment in line 2 of Algorithm 1, we obtain an affine control law that is always feasible 
and also optimal. 

When none of the trivial cases holds, we are in case [C4] of Section 4.1.2. Therefore, we can invoke the 
results from Lemma 4.3 to argue that the right side of the set Ar is exactly the set on line 7 of the algorithm, 
r-side(Ar) = {vq, ...,{),,}□ {vt} U {vr, . . . , ■Ufe}. In this setting, we can now formulate the first claim about 
system (37) and its solution: 

Lemma 4.4. System (37) is always feasible, and the solution satisfies: 

1- -b^<q^ <0, Vie {l,...,k}. 

2. L< q{w) <U, Vit) e Hk- 

Proof. Note first that system (37) has exactly fc + 3 unknowns, two for the cotangents Kjj, Kl, and one for 
each coefficient qi,0 < i < k. Also, since |r-side(Ar)| < |ext(Ar)| < k + 1, and there are exactly |r-side(Ar)| 
matching constraints, and fc + 3 — |r-sidc(Ar)| alignment constraints, it can be immediately seen that the 
system is always feasible. 

Consider any qi with i g {1, . . . , s} U {r + 1, . . . , fc}. From the matching conditions, we have that 
qi ~ u*{vi) — u*{vi_i). By property P3 from Section 2, the difference in the values of the optimal control 
law u*{-) satisfies: 

U*{V,) - U*{v,^i) = U*{92[V,]) - U*(^?2K-l]) 

(byP3) =-f-{92[v,]-92[v,^i]) 

''=-f-b,, where /e [0,1]. 

Since, by (32), bj > 0, Vj £ {1, ... , fc}, we immediately obtain —bi < < 0, for i g {1, . . . , s}U{r+l, . . . , fc}. 
Now consider any index i e {s + 1, . . . ,i Ar}, where tAr = min(t,r). From the conditions in system (37) 



15 



for alignment below t, we have qi 

t/\r 
i=s+l 



aj — Ku -bj 

Ku~c 



By summing up all such relations, we obtain: 



Etf\r \-\tf\T T 



<^ (using the matching) 



Kr. 



J2l=s+1 a-i+C - {u*{VtAr) - U*{Vs)) 





U*{VtAr) 


-[E-=o«' + c-^*(^.)] 






U*{VtAr) 





(29) 7l[ffArj - 7lFsj (?5) 



cotan({>s, {>tAr) ■ 



In the first step, we have used the fact that both Vs and Vmin{t.r) a-re matched, hence the intermediate 
coefficients qi must sum to exactly the difference of the values of at Vinin(t.r) ^md respectively. In 
this context, we can see that Ku is simply the cotangent of the angle formed by the segment [i's, 'Umin(t,r)] 
with the horizontal (i.e. 71) axis. In this case, we can immediately recall result (34) from Lemma 4.3, to 
argue that Kjj > Combining with (31) and (32), we obtain: 

a. + i (32) OniinCt r) CLt (31) 
^min(t,r) &t 

Therefore, we immediately have that for any i G {s + 1, . . . , min(i, r)}, 
\ Qi — Ku ■ bi <0 ai - Ku ■ b, 



Kr 



OO 



Ku - c 



<0 , 



c-bi > 
a,; 



c • b,. 



q^ + b^> 0. 



KjT - c 



The argument for indices i G {max(t, s) + 1, . . . ,?'} proceeds in exactly the same fashion, by recognizing 
that Ku defined in the algorithm is the same as cotan (■Umax(t,s)> Vr), and then applying (34) to argue that 



A L < ^ < ^ 



< 



at+l 



< c. This will allow us to use the same reasoning as above, completing the 



br — b„ 

proof of part {i) of the claim. 

To prove part (m), consider any w e TLk = [0,1]^ Using part (i), we obtain: 



q{w) = go 



^qt-Wi < (since e [0, l],gi < 0) < go u*{vo) < U , 

k 



qi-l = u*{vk) > L 



Note that in step (**), we have critically used the result from Lemma 4.3 that, when 8 ^ Buu, the points 
Vq, Vk are always among the points on the right side of Ar, and, therefore, we will always have the equations 
go = 'u*{vo),qo + EiLi 1i ~ ^*(''^fc) among the matching equations of system (37). For the last arguments, 
we have simply used the fact that the optimal control law, u*{-), is always feasible, hence L < u*{-) < U. □ 

This completes our first goal, namely proving that the affine controller q{w) is always robustly feasible. 
To complete the construction, we introduce the following final result: 



Lemma 4.5. Using the affine control law q(w) computed in Algorithm 1 satisfies the first equation in (35). 
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Proof. From (36), the affine controller q{w) induces the generators a + c ■ q and 6 + q for the zonogon F. 
This implies that F will be the Minkowski sum of the following segments in R^: 



ai + c - qi 




as +c- 






Ku ■ {bs+i + qs 


+ 1 


) ■ 






min{t,r] 


+ ^min(i,r)) 


bi+qi 


, . . . , 


bs + qs 




bs+i + qs+i 






, . . . , 




5min(i,r) 


■ (^max(t,s) + l 


^~ Qmax 


{t,s) + l) 




Kl ■ {br + qr) 






C • qr+1 




flfe + c ■ qk 


^max(t,s) + l + 


'Zmax(t,s) + 1 


, . . . , 


br + qr 


5 




br+l H 


- qr+1 


■)••'') 


bk + qk 



(38) 

From Lemma 4.4, we have that qi + bi > 0, Vi £ {1, . . . , fc}. Therefore, if we consider the points in R^: 

_ qj) I , Vi e {0,...,fc}, 

j=o 

we can make the following simple observations: 



Y,{aj+c-qj), ^(6j 

]=0 j=o 



• For any vertex Vi G Q, i G {0, . . . , k}, that is matched, i.e. Vi E r-side(Ar), if we let Wi represent the 
unique^ vertex of the hypercubc Hk projecting onto Vi, i.e. vi = {6i{wi),92{wi)), then we have: 

(36) , f . , (37) , , (29) . 

y. = (7i(u;,),72(u;j)) = {M'"i),l2{vi)) = v,. 

The first equality follows from the definition of the mapping that characterizes the zonogon F. The 
second equality follows from the fact that for any matched vertex Vi, the coordinates in F and F are 
exactly the same, and the last equality is simply the definition of the point Vi. 

• For any vertex vt G 0, j £ {0, . . . , fc}, that is not matched, we have: 

Vi e [Vs^ ymin(t,r)]. Vi £ {s + 1, . . . , min(t, r) - 1} 
Vi e [ymax(t.s)>yr]: £ {max(i, s) + 1, . . . ,r - 1}. 

This can be seen directly from (38), since the segments in R^ given by [y^, y^_|_i], . . . , [ymin(t.r)-ii ymin(t,r) 
are always aligned (with common cotangent, given by Ku), and, similarly, the segments [yijiax(t s) > ymax(t 
. . . , [Ur-i, Vr] s-re also aligned (with common cotangent K^). 

This exactly corresponds to the situation shown earlier in Figure 6. By combining the two observa- 
tions, it can be seen that the points ^y„, yi, . . . , J/^j ymax(t,s)) ymin(t,r)7 • • • , l/fc} wiU satisfy the following 
properties: 

yi=Vi, WviE r-side(Ar) , 
cotan (yp, y^) > cotan (yi, ya) > ' ' ' > cotan (y,_i, yj > cotan (y^, ymin(t,r)) > 
> cotan (y,„ax(t,s): Vr) > cotaii (y,,, y,,+i) > • • • > cotan (y^^i, y^) , 

where the second relation follows simply because the points Vi £ r-side(Ar) are extreme points on the right 
side of a convex hull, and thus satisfy the same string of inequalities. This immediately implies that this set 
of yj (corresponding to Vi £ r-sidc(Ar)) exactly represent the right side of the zonogon F, which, in turn, 
implies that F = z-huU ({vq, Vi, . . . , v^, Vj„ax(t,s), *mm(t,r), *r, I'r+i, • • ■ , t>fe}) • But then, by Corollary 4.1, 
the maximum value of problem (OPT) in (26) is equal to the maximum value of problem [AFF) in (35), 
and, since the former is always JmM, so is that latter. □ 

This concludes the construction of the affine control law q{w). We have shown that the policy computed 
by Algorithm 1 satisfies all the conditions in (35), i.e. is robustly feasible (by Lemma 4.4) and, when used in 
conjunction with the original convex state costs, preserves the overall optimal min-max value JmM (Lemma 
4.5). 



^This vertex is unique due to our standing Assumption 2 that the number of vertices in is 2k (also see part (iv) of Lemma 
7.2 in the Appendix). 
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4.3 Construction of the Afflne State Cost. 

Note that we have essentially completed the first part of the induction step. For the second part, we would 
still need to show how an afhne stage cost can be computed, such that constraints (13) and (14) are satisfied. 
We will return temporarily to the notation containing time indices, so as to put the current state of the 
proof into perspective. 

In solving problem {AFF) of (35), we have shown that there exists an affine qk+i{w''^^) such that: 
JraM = max [9i{w''+') + ck+i ■ qk+i{w''+^) + 9k+i {02{w''+') + qk+i{'^''+'))] 

^'i^ max [^,{w^+^)+gu+i{l2{w'^+'))]. 
Using the definition of .9fe+i(-) from (7), we can write the above (only retaining the second term) as: 



JraM = max 



7l(u;^+l) ^ [/lfc+2 (72 ('*«''+') + Wk+2) + Jfe+2(72('»'^+') + U.fc+2)] 

t«fc+l6Wfc + i 



max [ ^,{w'^+^) + {■K2{w^+^)) + JI+2 {^2{w^^^)) ] , 

where 7ri(it)'^+^) = 71(11;'"'+-'-), and 7r2(if'^+^) == 72(11''^+^) + Wk+2- Is is easy to note that: 

n^=*'(^i(«;^+2),^2(«^'+')) (39) 

represents yet another zonogon, obtained by projecting a hypcr-rectanglc TLk+i C R'"'+^ into R^. It has 
a particular shape relative to the zonogon T = (71,72), since the generators of 11 are simply obtained by 
appending a and a 1, respectively, to the generators of T , which implies that 11 is the convex hull of two 
translated copies of F, where the translation is occuring on the 112 axis. As it turns out, this fact will bear 
little importance for the discussion to follow, so we include it here only for completeness. 

In this context, the problem we would like to solve is to replace the convex function hk+2{'^2{w^^'^)) with 
an affinc function Zk+2{w^^'^) , such that the analogues of conditions (13) and (14) are obeyed: 

zfe+2(-u;^+') > hk+2{TT2{w''+^)), V-u;^+2 e Hk+i 

JraM = max [ T:i{w^+'^) + Zfc+2(u;'^'+^) + J^!+2(7r2(t«''+^)) ] . 

We can now switch back to the simplified notation, where all the time subscripts or superscripts are removed. 
Furthermore, to preserve as much of the familiar notation from Section 4.1.1, wc will denote the generators 
of zonogon tt by a, 6, so that wc have: 

fc+i fc+i 
7ri(u;) = ao + ^ fli • , t:2{w) ^ + ^^bi ■ Wi (40) 

i=l i=l 

In perfect analogy to our discussion in Section 4.1. we can introduce: 

= argmax < argmin {tt G 11} > ; = 20 — {O is the center of 11) (41) 

TTi I, 71-2 J 

Vq = , . . . , Vp^ , . . . , V2pi ~ (counter-clockwise numbering of the vertices of 11) . 

Without loss of generality, we will, again, work under Assumptions 1, 2, and 3, i.e. we will analyze the 
case when Hk+i = [0,l]'=+\ pi ^ k + 1 (the zonogon 11 has a maximal number of vertices), and Vi = 
[1, 1, . . . , 1, 0, . . . , 0] (ones in the first i positions). Furthermore, we will again use 7ri,2(if) for w G Ti-k+i to 
denote the mapping from Hk+i ^ 11 C R^, and 7ri.2[t;] to denote the coordinates of the point v G R^, and 
will use the shorthand notations h{vi), J* (vi) instead of h{TT2[vi]) and J*(7r2 [vj]), respectively. 
With the simplified notation, the goal is to find z{w) such that: 

z{w) > h{T:2{w)), Vw G Hk+i (42) 
max [tti + /i(7r2) + J*(7r2) ] = max [Tri{w) + z{w) + J* {Tr2{w))] (43) 
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In (43), the first maximization problem corresponds to the problem solved by the uncertainties, w, when 
the original convex state cost, h{'K2), is incurred. As such, the result of the maximization is always exactly 
equal to JmM, the overall min-max value. The second maximization corresponds to the problem solved by 
the uncertainties when the afiine cost, z{w), is incurred instead of the convex cost. Requiring that the two 
optimal values be equal thus amounts to preserving the overall min-max value, JmM, under the affine cost. 

Since h and J* are convex (see property P2 in Section 2), we can immediately use Lemma 4.1 to conclude 
that the optimal value in the first maximization problem in (43) is reached at one of the vertices Vq, . . . , Vk+i 
foimd in r-sidc(n). Therefore, by introducing the points: 

v,"^ {7:i[v,] + h{v,), ^2K]),Vze{0,...,fc + l}, (44) 

we can immediately conclude the following result; 

Lemma 4.6. The maximum in problem: 

{OPT) max [ + J* (^2) ] , 
where tti = tti + h{TT2), tt2 = 7^2, (tJ"!, 7r2) £ 11, 

is reached on the right side of: 

An = conv({'Uo,---,*fe+i})- (46) 



Proof. The result is analogous to Lemma 4.2, and the proof is a rehashing of similar ideas. In particular, 
first note that problem [OPT) is a rewriting of the first maximization problem in (43). Therefore, since the 
maximum of the latter problem is reached at the vertices Vi, i € {0, . . . ,k + 1}, oi zonogon 11, by using (44) 
we can conclude that the maximum in problem (OPT) must be reached on the set {vq, . . . ,Vk+i}. Noting 
that the function maximized in (OPT) is convex, this set of points can be replaced with its convex hull. An, 
without affecting the result [see Section 32 of Rockafellar, 1970]. Furthermore, since the function maximized 
is of the form ^1 + /(^2), with / convex, by applying the results in Corollary 4.1, and replacing the set by 
the right-side of its convex hull, r-side(An), the result of the optimization would remain unchanged. □ 



Continuing the analogy with the construction in Section 4.2, we rewrite the second optimization in (43) 

as: 

(AFF) max [ tti + J* (^2) ] , 

(*i,*2)6n ^^^^ 

where n=|(7ri,7r2) : tti{w) ""^ Tri{w) + z{w), T:2iw) tt2{w), w e Hk+i ^ ■ 

In order to examine the maximum in problem {AFF), we remark that its feasible set, 11 C R^, also represents 
a zonogon, with generators given by a + 2: and 6, respectively. Therefore, by Lemma 4.1, the maximum of 
problem {AFF) is reached at one of the vertices on the right side of II. 

Using the same key idea from the construction of the affine control law, we now argue that, if the 
coefficients of the affine cost, Zi, were computed in such a way that II represented the zonogon hull of the 
set of points {vq, . . . ,Vk+i}, then (by Corollary 4.1), the maximum value of problem (AFF) would be the 
same as the maximum value of problem {OPT). 

To this end, we introduce the following procedure for computing the affine cost z{w): 
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Algorithm 2 Compute afRnc stage cost z{w) 



Require: 'Ki{w),'K2{w), h{-), J* {■). 
1: Apply the mapping (44) to obtain Vi G {0, . . . , fc + 1}. 
2; Compute the set An = conv {{vq, . . . , Vk+i})- 

3; Let r-side(An) == {i's(i), • ■ • , ■Us(n)}, where s(l) < s(2) < ••• < s(n) G {0,...,fc + l} are the sorted 

indices of points on the right side of An • 
4; Solve the following system for Zj, (j G {0, . . . , fc + 1}), and -ftrs(i), (i G {2, . . . , n}): 

{zo + zi-\ h Zs(t) = h {vs(i)) , yvs{i) e r-side(An) (matching) 
= Vj G {s(*- + G {2,...,n}, (alignment) ^^^^ 

5: Return z(u;) = 20 + '^iH ' Wi- 



To visualize how the algorithm is working, an extended example is included in Figure 7. 



7r2 




Original zoiiotopc E. t^i Points Vj and j/^ e r-sidc(n). t^i 

Figure 7; Matching and alignment performed in Algorithm 2. 



The intuition behind the construction is the same as that presented in Section 4.2. In particular, the 
matching constraints in system (48) ensure that for any vertex w of the hypcrcubc Tik+i that corresponds 

to a potential maximizer in problem {OPT) (through w G Tik+i G 11 Vi G r-side(An)), the value 

of the affinc cost z{w) is equal to the value of the initial convex cost, h{vi) = h{TT2('w)), implying that the 
value in problem {AFF) of (47) at {Tti{w),TT2{w)) is equal to the value in problem (OPT) of (45) at Vi. 
The alignment constraints in system (48) ensure that any such matched points, (7ri(K;), 712(10)), actually 
correspond to the vertices on the right side of the zonogon II, r-side(II), which implies that, as desired, 
n = z-huU ({{jo, . . . , Vk+i}). 

We conclude our preliminary remarks by noting that system (48) does not directly impose the robust 
domination constraint (42). However, as we will soon see, this result is a byproduct of the way the matching 
and alignment are performed in Algorithm 2. 



4.3.1 Affine Cost z( ) Dominates Convex Cost /i( ) and Preserves Overall Objective. 

In this section, we will prove that the afline cost z(w) computed in Algorithm 2 not only robustly dominates 
the original convex cost (42), but also preserves the overall min-max value (43). 
The following lemma summarizes the first main result: 

Lemma 4.7. System (48) is always feasible, and the solution z(w) always satisfies equation (43). 
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Proof. We first note that s(l) = and s{n) = k + 1, i.e. Vo,Vk+i G r-side(An). To see why that is 
the case, note that, by (41), Vq will always have the smallest 7r2 coordinate in the zonogon H. Since 
the transformation (44) yielding Vi leaves the second coordinate unchanged, it will always be true that 
Vq ~ a-rgmax^^ argmin^^ {vi, i S {0, . . . ,k + 1}}, which will immediately imply that Vq G r-side(An). The 
proof for Vk+i follows in an identical matter, since Vk+i has the largest 7r2 coordinate in 11. 
It can then be checked that the following choice of Zi always satisfies system (48): 



zq = h{vo); Zj 



s{i) 



\/j e {s{t-l) + l,...,s{i)}, Vz e {2,...,n}, 



^s(i-i)+i + • • • + + a,j(i_i)_|_i + • • • + as(j) ^ h{vs{i)) - h{vs[i-i)) + a<;(j_i)+i + • • • + a^fj) 



^s('i-l) + l + • • • + 6s(j) 



The proof of the second part of the lemma will be analogous to that of Lemma 4.5. To start, consider 
the feasible set of problem [AFF] in (47), namely the zonogon 11, and note that, from (40), its generators 
are given by a + z, b: 



a + z 




ai + 


b 




bi 



as{i) + Zs(i) a5(i)+i + Zs(i)+i . . . ttk+i + Zk+i 



By introducing the following points in R : 



(49) 



we have the following simple claims: 



• For any Vi G 11 that is matched, i.e. Vi G r-side(An), with Wi = [1,1,... 
unique^ vertex oiTik+i satisfying {TTi{wi),Tr2{wi)) = Vi, we have: 

Vl * = (TTliWi) + z{Wi),Tr2iWi)) ^^'^ (iTllVi] + h{Vi),TT2[Vi]) 



, 1, 0, . . . , 0] denoting the 



(44) , 
= Vi 



The first equality follows from the definition of the zonogon ft, the second follows because any Vi G 
r-side(An) is matched in system (48), and the third equality represents the definition of the points Vi. 

• For any vertex vj G 11, that is not matched, i.e. vj ^ r-side(An), and s(i) < j < s{i + 1) for some i, we 
have Uj G [ys(i) , ys{i+i)]- This can be seen by using the alignment conditions in system (48) in conjunc- 
tion with (49), since the segments inR^ given by [ys{z),ys{t)+i], [y.s(»)+i> y.s(<)+2]> 
will always be parallel, with common cotangent given by Kg(_i^iy 

For a geometric interpretation, the reader is referred back to Figure 7. Corroborating these results with the 
fact that {■Us(i), . . . , ■Ws(n)} = r-sidc(An) must always satisfy: 

cotan (■Us(i), Vs{2)) > cotan {vs(2), Vs{3)) > ■ ■ ■> cotan Vs(n)) , (50) 

we immediately obtain that the points 2/3(2); ■ • • iys(n)| (corresponding to Vi G r-side(An)) exactly 

represent the right side of the zonogon ft, which, in turn, implies that ft = z-huU {{vq, Vk+i})- But 

then, by Corollary 4.1, the maximum value of problem (OPT) in (45) is equal to the maximum value of 
problem (AFF) in (47), and, since the former is always JmM, so is that latter. □ 

In order to complete the second step of the induction, we must only show that the robust domination 
constraint (42) is also obeyed: 

z{w) > h (7r2(u;)) ^ zo + zi • wi H h Zfc+i • Wfe+i > h (bo + bi ■ wi ^ h bk+i ■ Wk+i) , Vii; G Hk+i- 

The following lemma will take us very close to the desired result: 



''Recall that we are, again, working under Assumption 2, which implies uniqueness by part (iv) of Lemma 7.2 in the Appendix) 
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Lemma 4.8. The coefficients for the affine cost z{w) computed in Algorithm 2 will always satisfy the following 
property: 



h (bo + fej(i) H h &j(m)) < ^0 + 2^(1) H K Zj(,„), VjXI), • . -.jim) £ {1, . . . , fc + l},Vm G {1, 



!}• 



Proof. Before proceeding with the proof, we will first list several properties related to the construction of the 
affine cost. We claim that, upon termination, Algorithm 2 will produce a solution to the following system: 



zo + zi + 



i zo + zi 

zi+ai 
bi 



Zs(2) 



h («s(n)) 



f'3(i.-l) + l 



Z3(2)+Q3(2) 



(51) 



K 



s{2} 



s(n) 



h{Vj)~h{Vo)+ai+---+aj , ^ , /i(t!3(2) ) CUj )+aj + l + ■ ■ ■+03(1) 
6i + -+b, ^ ^.^2) < 



fcj + l + ---+b3 



VjG{l,...,5(2)-l} 



/t(t>j)-/t(t>3(n-l))+a3(„^l) + l+---+aj 
6s(„_i) + iH h&j 



(52) 



(53) 



< K 



s(n) 



< 



h(D3(„))-ft(l'i)+ai + lH ha,(„) 

fcj + l + ---+fc3(„) ■ 



Vj € {s(ri- 1) + 1,. ..,s(n) - 1} 



Let us explain the significance of all the equations. (51) is simply a rewriting of the original system (48), 
which states that at any vertex fsji), the value of the affine function should exactly match the value assigned 
by the convex function /i(-), and the coefficients Zi between any two matched vertices should be such that 
the resulting segments, [zj + aj,6j], are aligned (i.e. the angle they form with the tti axis has the same 
cotangent, specified by Ki^.) variables). We note that we have explicitly used the fact that s(l) = 0, which 
we have shown in the first paragraph of the proof of Lemma 4.7. 

Equation (52) is a simple restatement of (50), that the cotangents on the right side of a convex hull must 
be decreasing. 

Equation (53) is a direct consequence of the fact that {^8(1)1 i's(2)j • ■ • i''^s(n)} represent r-side(An). To 
see why that is, consider an arbitrary j G {s{i) + 1, . . . , s(i + 1) — 1}. Since Vj ^ r-side(An), we have: 



~\ N (40), (44 

cot'An[ys(i), Vj) < cotan(t;j, ^ 

s{i)+i H h flj + jvj) - h {vsji)) ^ gj+i + ^ a,s(*+i) + h (v 



as(t}+i H ^Oj + h (Vj) - h (t's(j)) 



< K 



six 



< 



bj+1 + 



(»+i); 



h V 



s(i+l)J 



s(i+l) 



where, in the last step, we have used the mediant inequality® and the fact that, from (51), Ks{2) = 
cotan (^sfj), Vs(i+i)) = h(v,^,j) ^^^^^^ back to Figure 7 for a geometrical in- 

terpretation). 

With these observations, we will now prove the claim of the lemma. The strategy of the proof will be 
to use induction on the size of the subsets, m. First, we will show the property for any subset of indices 
j(l), . . . ,j{m) G {s(l) = 0, . . . , s(2)}, and will then extend it to j(l), . . . , j(™) G {«(«) + 1, . . . , s(i + 1)} for 
any i, and then to any subset of {1, . . . ,k + 1}. 



^If b,d>0 and^ < then i < < ^. 
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The following implications of the conditions (51), (52) and (53), are stated here for convenience, since 
they will be used throughout the rest of the proof: 

h{vs(^i^) ^ h{vo) ^ zn; /i(t;^(2)) = zo + zi H \~ Zs(^2)- (54) 

h{vj)-h{vo) <zi + --- + zj, VjG{l,...,s(2)-l}. (55) 

fl<...<fl<...<ff(2), Vje{l,...,.(2)-l}. (56) 
bi Oj bs(2) 



Their proofs are straightforward. (54) follows directly from system (51), and: 

h{vj) - h{va) + ai H h (^J) (g^) ^ h z,- + ai H ha,- 

< Ks(2) = ^ l--- + hj ^(55) true. 

(56) true. 



+zi _ . . . _ a-i+Zj _ . . . ^ °°(2)- 



'(51)- ai+zi — — aj+Zj _ ^ as(2)+Zs(2} 

n zonogon ^ > ■ ■ ■ > t- > ■ ■ ■ > 



We can now proceed with the proof, by checking the induction for m = 1. We would like to show that: 

h{bo + b,) <zo + Zj, e {l,...,s(2)} 



Writing bo + bj as: 



6o + 6j = (1 - A) • 6o + A ■ (6o + ■ • ■ + bj) 



bi + --- + bj 
we obtain: 

h{bo + bj) < (1 - A) • h{bo) + A • h{ba + --- + bj) 

= h(V,) 

= h{vQ) + ^-^ [h{vj) - h{vn) ] < (by (54) if j ^ s{2) or (55) otherwise) 

Oi + • • • + Oj 

< zq + ^ — (zi + • ■ • + Zj) < (by (56) and the mediant inequality) 

Oi + ■ ■ • + Oj 

< zo + Zj. 

Assume the property is true for any subsets of size m. Consider a subset j(l), . . . , j(m), j(m + 1), and, 
without loss of generality, let j(m + 1) be the largest index. With the convex combination: 

&* = feo + bj^i) + ■ • • + fej(,„) + &j(,n+l) 

= (1 - A) • (6o + bj(if H h 6j(„i)) + \- {bo + bi-\ h 6j(,„+i)-i + 6j(,„+i)), 



where A = 



b 



'i(m+l) 



(6l + 62 H h &j(m+l)) - (&i(l) + fej(2) H (- &j(m)) 

we obtain: 

h{b*) < (1 — A) • h{bo + fej(i) + ■ • • + fej(m)) + A ■ /i (^^^(m+i)) < (by induction hypothesis and (54), (55)) 
< (1 - A) ■ {zq + Zj(i) H h Zj(„i)) + X ■ {zq + zi -\ h 2;i(,„+i)) 



2^0 + Zj(l) H K Zji^jn) + 



bj{7n+l) 



(61 + 62 H \- bj(„i+i)) ~ (&j(i) + &j(2) H h 6j(m)) 

• [(zi + Z2 H h Zj(„;+i)) - (zj(i) + Zj(2) H 1- Zj(„))] < (by (56) and mediant inequality) 

< zo + Zj(i) H h Zj(™) + Zj(,„_,_i). 
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We claim that the exact same procedure can be repeated for a subset of indices from {s{i) + 1, . . . , s(i + 1)}, 
for any index i G {1, . . . ,n — 1}. We would simply be using the adequate inequality from (53), and the 
statements equivalent to (54), (55) and (56). The following results would be immediate: 

h {{bo + bi-\ h fes(,)) + H h < {za + zi + h + H + zj^^^-j, (57) 

Vie {l,...,7i}, Vj(l),...,j("^) e {s(^) + l,...,s(^ + l)}. 

Note that instead of the term bo for the argument of h{-), we would use the complete sum &o + + ■ ■ ■ + ^s(i) j 
and, similarly, instead of zq wc would have the complete sum zq + + • ■ ■ + Zs(i)- With these results, we 
can make use of the increasing increments property of convex functions: 

h{xi + A) - h{xi) ^ h{x2 + A) - h{x2) , 

< : , VA > U,Xi < X2 , 



A 



A 



to obtain the following result: 



bid 



(1) 



V 



-h{bo + + • • • + 



< 



< h 



j(-)6{l...,s(2)} i(-)6{s(2) + l,...,s(3)}/ 

/ 

bo+ bi-\ h5,(2) +&i(i)H 1- &i(/) I - ^ (&o + ^1 H ^bs(2)) < 

V 



(54), (57) 



all indices in {l,...,s(2)} 

< (zo + 21 H K Zs(2)) + 2i(l) H 'r Zi(^i) - (zo + Zl + 

= Z,(l) H h Zj(;) ^ 



=''''(l'a(2)) 



2s(2) 



(57) 



h {bo + &j(i) H h &j(,„) + 5i(i) + h < /i (&o + &i(i) H h + Zi(i) H h z^^ < 

< Zo + Zj(i) + h Zj(„) + Zj(i) H + Zi(;). 

We showed the property for indices drawn only from the first two intervals, {s(l) + 1, . . . , s(2)} and {s(2) + 
l,...,s(3)}, but it should be clear how the argument can be immediately extended to any collection of 
indices, drawn from any intervals. We omit the details for brevity, and conclude that the claim of the lemma 
is true. □ 

We are now ready for the last major result: 

Lemma 4.9. The affine cost z{w) computed by Algorithm 2 always dominates the convex cost h{'K2{w)): 

/ fe+1 \ fc+i 



h\bo + '^b.i-w, \ < zq + ^ Zi -Wi, Vu; e Hk+i = [0, 1] 



fe+1 



Proof. Note first that the function f{w) h ^6o + 'I2i=i bi ■ Wij — (zq + X]i=^i ■^i ' '^i) ^ convex function 
of w. Furthermore, the result of Lemma 4.8 can be immediately rewritten as: 



fe+1 



fe+1 



/i Uo + ^6, -lyj < zo + ^z, -u;,, Vio G {0,l}'=+i ^ f{w) <0,Wwe {0,1} 



fe+1 



1=1 



Since the maximum of a convex function on a polytope occurs on the extreme points of the polytopc, and 
ext(?ifc+i) = {0,1}*^+^, wc immediately have that: maxweHk+i f{''^) = ^'^^^we{o.i}''+^ f{'^) ^ 0, which 
completes the proof of the lemma. □ 

Wc can now conclude the proof of correctness in the construction of the affine stage cost, z{w). With 
Lemma 4.9, wc have that the affine cost always dominates the convex cost /i(-), thus condition (42) is obeyed. 
Furthermore, from Lemma 4.7, the overall min-max cost remains unchanged even when incurring the affine 
stage cost, z{w), hence condition (43) is also true. This completes the construction of the affine cost, and 
hence also the full step of the induction hypothesis. 
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4.3.2 Proof of Main Theorem. 



To finalize the current section, we summarize the steps that have lead us to the result, thereby proving the 
main Theorem 3.1. 

Proof of Theorem 3.1. In Section 4.1, we have verified the induction hypothesis at time k = \. With the 
induction hypothesis assumed true for times t = 1, . . . , fc, we have listed the initial consequences in Lemma 
4.1 and Corollary 4.1 of Section 4.1.1. By exploring the structure of the optimal control law, 
and the optimal value function, J^^j^(xfc+i), in Section 4.1.2, we have finalized the analysis of the induction 
hypothesis, and summarized our findings in Lemmas 4.2 and 4.3. 

Section 4.2 then introduced the main construction of the affine control law, qk+i{w^^^), which was 
shown to be robustly feasible (Lemma 4.4). Furthermore, in Lemma 4.5, we have shown that, when used in 
conjuction with the original convex state costs, ft.fe+i {xk+2)i this affine control preserves the min-max value 
of the overall problem. 

In Section 4.3, we have also introduced an affine stage cost, Zk+2{w^^'^): which, if incurred at time fc + 1, 
will always preserve the overall min-max value (Lemma 4.7), despite being always larger than the original 
convex cost, hk+i {xk+2) (Lemma 4.9). □ 

4.3.3 Counterexamples for potential extensions. 

On first sight, one might be tempted to believe that the results in Theorem 3.1 could be immediately extended 
to more general problems. In particular, one could be tempted to ask one of the following natural questions: 

1. Would both results of Theorem 3.1 (i.e. existence of affine control laws and existence of affine stage 
costs) hold for a problem which also included linear constraints coupling the controls ut across different 
time-steps? (see Bcn-Tal et al. [2005b] for a situation when this might be of interest) 

2. Would both results of Theorem 3.1 hold for multi-dimensional linear systems? (i.e. problems where 
Xk e R'', Vfc, with d>2) 

3. Are affine policies in the disturbances optimal for the two problems above? 

In the rest of the current section, we would like to show how these questions can all be answered negatively, 
using the following simple counterexample: 

7" = 4, Cfc = 1, hk{xk+i) = max{18.5 ■ x^+i, -24 • Xk+i], = 0, Uk = 00, Vfc G {1, . . . , 4} , 
{CEx) wi e [-7,0],W2 e [-11,0],W3 e [-8,0],W4 e [-44,0] , 

k 

< 10- fc ,Vfc G {1,...,4}. 

i=l 

The first two rows describe a one-dimensional problem that fits the conditions of Problem 1.1 in Section 1. 
The third row corresponds to a coupling constraint for controls at different times, so that the problem fits 
question (i) above. Furthermore, since the state in such a problem consists of two variables (one for Xk and 
one for X^iLi "^k), the example also fits question (ii) above. 

The optimal min-max value for the counterexample (CEx) above can be found by solving a stochastic 
optimization problem (sec Ben-Tal et al. [2005b]), in which non-anticipatory decisions are computed at all 
the extreme points of the uncertainty set, i.e. for {wi ,Wi} x ,^2} x {w^ ,^3} x {w^ , W4}. The resulting 
model, which is a large linear program, can be solved to optimality, resulting in a corresponding value of 
approximately 838.493 for problem {CEx). 

To compute the optimal min-max objective obtained by using affine policies qkiw'^) and incurring affine 
costs Zk{w''^^), one can ammend the model {AARC) from Section 3 by including constraints for the cumu- 
lative controls (see Ben-Tal et al. [2005b] for details), and then using (15) to rewrite the resulting model as 
a linear program. The optimal value of this program for counterexample (CEx) was approximately 876.057, 
resulting in a gap of 4.4%, and thus providing a negative answer to questions (i) and (ii). 

To investigate question (iii), wc remark that the smallest objective achievable by using affine policies 
of the type qk{w^) can be found by solving another stochastic optimization problem, having as decision 
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variables the afBne coefficients {qk,t}o<t<k<T, as well as (non-anticipatory) stage cost variables z^f for every 
time step k G {1, • ■ • ,T} and every extreme point w of the uncertainty set. Solving the resulting linear 
program for instance {CEx) gave an optimal value of 873.248, so strictly larger than the (true) optimum 
(838.493), and strictly smaller than the optimal value of the model utilizing both afhnc control policies and 
affine stage costs (876.057). 

Thus, with question (iii) also answered negatively, we conclude that policies that arc afhnc in the dis- 
turbances, qk(w^), are in general suboptimal for problems with cumulative control constraints or multiple 
dimensions, and that replacing the convex state costs hk{xk+i) by (larger) affine costs Zk{w''~^^) would, in 
general, result in even further deterioration of the objective. 

5 An application in inventory management. 

In this section, we would like to explore our results in connection with a classical inventory problem. This 
idea was originally introduced by Ben-Tal et al. [2005b], in the context of a more general model: the retailer- 
supplier with flexible commitment contracts problem (RSFC). We will first describe a simplified version of 
the problem, and then draw a very interesting connection with our results. 

The setting is the following: consider a single-product, single-echelon, multi-period supply chain, in 
which inventories are managed periodically over a planning horizon of T periods. The unknown demands 
Wt from customers arrive at the (unique) echelon, henceforth referred to as the retailer, and are satisfied 
from the on-hand inventory, denoted by Xt at the beginning of period t. The retailer can replenish the 
inventory by placing orders Ut, at the beginning of each period t, for a cost of Ct per unit of product. These 
orders are immediately available, i.e. there is no lead-time in the system, but there are capacities on how 
much the retailer can order: Lt < Ut < Ut- After the demand Wt is realized, the retailer incurs holding 
costs Ht ■ max{0,xt + ut — wt\ for all the amounts of supply stored on her premises, as well as penalties 
Bt • max{wt — xt — ut, 0}, for any demand that is backlogged. 

In the spirit of robust optimization, we will assume that the only information available about the demand 
at time t is that it resides within a certain inverval centered around a nominal (or mean) demand dt , which 
results in the uncertainty set Wt = { |ift — (it | ^ P' dt}, where p S [0,1] can be interpreted as an uncertainty 
level. As such, if we take the objective function to be minimized as the cost resulting in the worst-case 
scenario, we immediately obtain an instance of our original Problem 1.1, with at = /3t = 1, 7f = — 1, and the 
convex state costs ht{-) denoting the Newsvendor cost, ht{xt+i) = Ht ■ maxjxt + ut — Wt,0} + Bt ■ max{wt — 

Xt - Mt,0}. 

Therefore, the results in Theorem 3.1 are immediately applicable to conclude that no loss of optimality 
is incurred when we restrict attention to afRne order quantities qt that depend on the history of available 
demands at time t, qt{w^) = fZt,o + X]t=i 1^,^ ' ""^t; and when we replace the Newsvendor costs ht{xt+i) by 
some (potentially larger) affine costs zt{w*^^). The main advantage is that, with these substitutions, the 
problem of finding the optimal afhne policies becomes an LP (see the discussion in Section 3 and the paper 
by Ben-Tal et al. [2005b] for more details). 

The more interesting connection with our results comes if we recall the construction in Algorithm 1. In 
particular, we have the following simple claim: 

Proposition 5.1. // the affine orders qt{w^) computed in Algorithm 1 are implemented at every time step 

t, and we let: Xk{w^) = xi + X]t=i ilti'^^) ~ Wt) = Xtfl + X)t=i •'^fe,* ' denote the affine dependency of 
the inventory Xk on the history of demands, , then: 

1. If a certain demand Wt is fully satisfied by time k > t + 1, i.e. Xk,t = 0, then all the (affine) orders q^ 
placed after time k will not depend on Wt . 

2. Every demandwt is at most satisfied by the future orders qk, k > t+l, and the coefficient qk^t represents 
what fraction of the demand Wt is satisfied by the order qk ■ 

Proof. To prove the first claim, recall that, in our notation from Section 4.1.1, Xk ^ 62 = bo + X]t=i ' "^t- 
Applying part (i) of Lemma 4.4 in the current setting^, we have that < qk.t < ^Xk^t- Therefore, if 

®The signs of the inequalities are changed because every disturbance, wt, is entering the system dynamics with a coefBcicnt 
— 1, instead of as was the case in the discussion from Section 4.1.1. 
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Xk t = 0, then qk t = 0, which imphcs that Xk+i t = 0. By induction, wc immediately get that t = 0, Vt G 
{k,...,T). 

To prove the second part, note that any given demand, wt, initiaUy has an affine coefficient of —1 in the 
state Xt+i, i.e. Xt+i,t = ^ 1- By part [i) of Lemma 4.4, < qt+i,t < —Xt+i,t = 1, so that qt+i,t represents a 
fraction of the demand Wt satisfied by the order qt+i- Furthermore, Xt+2,t = Xt+i^t + <lt+i,t £ 0]i so, by 
induction, wc immediately have that q^^t & [0, 1], Vfc > t + 1, and X]fc=t+i Ik.t < 1. CH 

In view of this result, if wc think of {qk}k>t+i as future orders that are partially satisfying the demand 
Wt, then every future order quantity qk{'w^) will satisfy exactly a fraction of the demand Wt (since the 
coefficient for wt in will always be in [0, 1]), and every demand will be at most satisfied by the sequence 
of orders following after it appears. This interpretation bears some similarity with the unit decomposition 
approach of Muharremoglu and Tsitsiklis [2007] , where every unit of supply can be interpreted as satisfying 
a particular unit of the demand. Here, we are accounting for fractions of the total demand, as being satisfied 
by future order quantities. 



6 Conclusions. Future Directions. 

We have presented a novel approach for theoretically handling robust, multi-stage decision problems. The 
method strongly utilized the connections between the geometrical properties of the feasible sets (zonogons), 
and the objective functions being optimized, in order to prune the set of relevant points and derive properties 
about the optimal policies for the problem. We have also shown an interesting implication of our theoretical 
results in the context of a classical problem in inventory management. 

On a theoretical level, one immediate direction of future research would be to study systems with mixed 
(polyhedral) constraints, on both state and control at time t. Furthermore, we would like to explore the 
possibility of utilizing the same proof technique in the context of multi-dimensional problems, as well as for 
more complicated uncertainty sets W. 

Second, we would like to better understand the connections between the matching performed in Algorithm 
2 and the properties of convex (and supermodular) functions, and explore extensions of the approach to 
handle cost functions that are not necessarily convex, as well as non- linear cost structures for the control ut. 
Another potential area of interest would be to use our analysis tools to quantify the performance of affine 
policies even in problems where they arc known to be suboptimal (such as the one suggested in Section 
4.3.3). This could potentially lead to fast approximation algorithms, with solid theoretical foundations. 

On a practical level, we would like to explore potential applications arising in robust portfolio optimiza- 
tion, as well as operations management. Also, we would like to construct a procedure that mimics the 
behavior of our algorithms, but does not require knowledge of the optimal value functions J*(-) or optimal 
controllers u*{-). One potential idea would be to explore which types of cuts could be added to the linear 
program {AARC), to ensure that it computes a solution as "close" to the affine controller q{w) as possible. 



7 Appendix. 

7.1 Dynamic Programming Solution. 

This section contains a detailed proof for the solution of the Dynamic Programming formulation, initially 
introduced in Section 2. Recall that the problem wc would like to solve is the following: 



J dot 

JmM = mni 



Ci ■ ui + max 

wi 



hi{x2) + ' 



mm 



Ck ■ Uk + max [hk{xk+i) + 



(DP) 



-niin ct • ut + max /iT(^T+i) 



s.t. Xk+1 = Xk + Uk + Wk 
Lk <Uk < Uk 
Wk eWk^ [wk ,Wk] 



Vfce {i,2,...,r}. 
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which gives rise to the corresponding BeUman recursion: 



Jki^k) min 

Lk<Uk<Uk 



Ck-Uk+ max [ hk{xk + Uk + Wk) + Jl+i {xk + Uk + Wk) 



According to our definition of running cost and cost-to-go, tire cost at T + 1 is J^+i = 0, which yields the 
following Bellman recursion at time T: 



J^{xt) min 



ct ■ ut + max hrixT + ut + wt) 



First consider the inner (maximization) problem. Letting yx == xt + ut, we obtain: 

gTivx) = max hxixT + ut + wt) 

(since /it(') convex) = max{/ir [yx + w^.) , hx (j/t + wt)} (58) 

Note that gr is the maximum of two convex, coercive functions of j/t, hence it is also convex and coercive 
(see Theorem 5.5 in RockafcUar [1970] or Proposition 1.2.4 in Bertsekas [2003] for details). The outer 
(minimization) problem at time T becomes: 

J^{xt) ~ min • ut + grixr + ut) 

= —Ct ■ XT + min [ ct ■ {xt + ut) + gTixT + ut) ] 

For any xt, ct • {xt + ut) + gT{xT + "t) is a convex function of its argument yT = xt + ut- As such, by 
defining [y^ ,yf] to be the compact set of minimizers of the convex and coercive function ct ■ y + griy), we 
obtain that the optimal controller and optimal value function at time T will be: 

Ut, if XT <y_r^- Ut 

Ut{xt) — { —XT + y*, otherwise (59) 
Lt, if Xt > Vt ~ Lt 

Ct -Ut + gT{xT + Ut), if xt < y^ - Ut 

Jt{xt) = { Ct ■ {y* ~ xt) + gTiy*), otherwise (60) 

^ct ■ Lt + gT{xT + Lt), if xt > Vt ~ Lt 
where y* e [y^ ^Ijt] ■ 

The following properties are immediately obvious: 

1. u^{xt) is piecewisc affine (with at most 3 pieces), continuous, monotonically decreasing in xt- 

2. J^{xt) is convex, since it represents a partial minimization of a convex function with respect to one 
of the variables (see Proposition 2.3.6 in Bertsekas [2003]). 

The results can be immediately extended by induction on k: 

Lemma 7.1. The optimal control policy M^(xfc) is piecewise affine, with at most 3 pieces, continuous, and 
monotonically decreasing in Xk- The optimal objective function J^ixk) is convex in Xf. 

Proof. The induction is checked at fc = T. Assume the property is true at fc + 1. Letting yk Xk + Uk, the 
Bellman recursion at k becomes: 

Jki^k)= min [ck-Uk+gkixk+Uk) ] 

Lk<Uk<Uk 

9k iVk) max [ hk{yk + Wk) + Jk+i (Vk + Wk) ] . 
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Consider first the maximization problem. Since hk is convex, and (by the induction hypothesis) J^+i is 
also convex, the maximum will be reached on the boundary of Wfe — [w/. ,Wk]'- 

9k {Vk) = max { hk(yk + w^) + Jfc+i ivk + Wk) : hk{yk + Wk) + J^+i {yu + Wk) } (61) 

and gkivk) will be also be convex. The minimization problem becomes: 

■Jki^k) ^ min [ck-Uk+gk{xk + Uk) ] 

Lk<Uk<Uk 

= -'Ck-Xk+ min [ Ck ■ (xk + Uk) + gk {xk + Uk) ] (62) 

Lfc<Ufc<C/fc 

Defining, as before, [y , y^.] as the set of minimizers of Ck ■ y + gk{y), we get: 







if Xfe < 


-Uk 




1 -a;fe + 


?/*, otherwise 






[Lk. 


if Xfe > Ijk 


-Lk 




fcfe ■ C/fe 


+ gk{xk + Uk), 


if Xk <y_^- Uk 




1 Cfe • (y* 


-xk) + gk{y*), 


otherwise 




[cfe • Lk 


+ gk{xk + Lk), 


if a;fe > Ijk - Lk 



where y* G [y^ , y^,] 



(63) 



(64) 



In particular, u^. will be piecewise afhne with 3 pieces, continuous, monotonically decreasing, and will 
be convex (as the partial minimization of a convex function with respect to one of the variables). A typical 
example of the optimal control law and the optimal value function is shown in Figure 1 of Section 2. □ 



7.2 Zonotopes and Zonogons. 

In this section of the Appendix, we would like to outline several useful properties of the main geometrical 
objects of interest in our exposition, namely zonotopes. The presentation here parallels that in Chapter 7 of 
Ziegler [2003], to which the interested reader is referred to for a much more comprehensive treatment. 

Zonotopes are special polytopes that can be viewed in various ways: as projections of cubes, as Minkowski 
sums of line segments, and as sets of bounded linear combinations of vector configurations. Each description 
gives a different insight into the combinatorics of zonotopes, and there exist some very interesting results that 
unify the different descriptions under a common theory. For our purposes, it will be sufficient to understand 
zonotopes under the first two descriptions. In particular, letting Tlk denote the fc-dimensional hypercube, 
Hk = {w en'' : <w, < 1, Vi}, we can introduce the following definition: 

Definition 7.1 (7.13 in Ziegler [2003]). A zonotope is the image of a cube under an affine projection, that 
is, a d-polytope Z C R"^ of the form 

Z = Z(V):=V -Hk + z^iVw + z : w eUk} 

k 

= {a;eR'^ : x ^ z + WiVi, < < 1} 

i=l 

for some matrix (vector configuration) V — {vi, . . . ,Vk) G R''^'^'. 

The rows of the matrix V are often referred to as the generators defining the zonotope. An equivalent 
description of the zonotope can be obtained by recalling that every fc-cube Tik is a product of line segments 
T^k = Ti-i X ■ ■ ■ xTCi- Since for a linear operator tt we always have: 7r(7ii x • • • x Hi) = tt{Hi) + • • • + 7r(?ii), 
by considering an affine map given by tt{w) = Vw + z, it is easy to see that every zonotope is the Minkowski 
sum of a set of line segments: 

Z{V)^[0,Vi] + --- + [0,Vp]+z. 
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For completeness, we remark that there is no loss of generality in regarding a zonotope as a projection 
from the unit hypercube Tife, since any projection from an arbitrary hyper-rectangle in R*"' can always be 
seen as a projection from the unit hypercube in R*^'. To see this, consider an arbitrary hyper-rectangle in 

Wk = [Wi , Wi] X , W2] X • • • X [w,. , Wfe], 

and note that, with V G R"^^*^, and £ R'' denoting the j-th row of V, the j-th component of Z{V) = V ■ 
Wk + z can be written: 



z(y)j = Zj + ^ (oj • wi) 



^ Oi • + X] ' ~ — i) ■ "^here G [0, 1], V 1 < i < /c. 



An example of a subclass of zonotopes are the zonogons, which are all centrally symmetric, 2-dimensional 
2p-gons, arising as the projection of p-cubes to the plane. An example is shown in Figure 2 of Section 4.1. 
These are the main objects of interest in our treatment, and the following lemma summarizes their most 
important properties: 



Lemma 7.2. Let Tik = [0,1]^^ be a k- dimensional hypercube, k > 2. For fixed a,b <E R*"' and ao,6o G R, 

r «T 

consider the affine transformation tt : R''' R^, Tr{w) 



a 



an 
bo 



and the zonogon Q C R 



e = 7r(Wfe) {0 G R2 : 3t(; G Uk s.t. ^ t:(w)] . 
If we let Ve denote the set of vertices of 0, then the following properties are true: 

1. 3 O G such that Q is symmetric around O : Va; G => 20 — cc G O. 

2. |Ve| = 2p < 2k vertices. Also, p < k if and only if3i=^jE {1, . . . , fc} such that rank 
2. 



hi hi 



< 



3. If we number the vertices of Ve in cyclic order: 

Ve ^ [VQt ■ ■ ,Vi,Vi+i, . . . ,V2p-l) iv2p+i V(^2p+i} mod (2p)) 

Vi^p, and we have the following representation for Q as a Minkowski sum of line 



then 20 — Vi 
segments: 



e = o 



2 ' 2 

Vi - Vj-l 



-1 < A, < 1. 



4- If 3wi,W2 G Ti-k such that Vi = 7t{wi) = V2 == t:{w2) and Vi_2 G Ve, then 3j G {1, . . .,k} such that 
«j = = 0. 

5. With the same numbering from (Hi) and k — p, for any i G {0,...,2p — 1}, the vertices of the 
hypercube that are projecting to Vi and fi+i, respectively, are adjacent, i.e. they only differ in exactly 
one component. 

Proof. We will omit a complete proof of the lemma, and will instead simply suggest the main ideas needed 
for checking the validity of the statements. 

For part (i), it is easy to argue that the center of the hypercube. On = [1/2, 1/2, . . . , 1/2]-'", will always 
project into the center of the zonogon, i.e. O — t:{Ou)- This implies that any zonogon will be centrally 
symmetric, and will therefore have an even number of vertices. 
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Part (m) can be shown by induction on the dimension k of the hypercube, Hk- For instance, to prove 
the first claim, note that the projection of a polytope is simply the convex hull of the projections of the 
vertices, and therefore projecting a hypercube of dimension k simply amounts to projecting two hypcrcubcs 
of dimension fc — 1, one for = and another for Wk — I, and then taking the convex hull of the two 
resulting polytopes. It is easy to see that these two polytopes in are themselves zonogons, and are 
translated copies of each other (by an amount [ak,bk]'^). Therefore, by the induction hypothesis, they have 
at most 2(fc — 1) vertices, and taking their convex hull introduces at most two new vertices, for a total of at 
most 2{k — 1) + 2 = 2fc vertices. The second claim can be proved in a similar fashion. 

One way to prove part (iii) is also by induction on p, by taking any pair of opposite (i.e. parallel, of the 
same length) edges and showing that they correspond to a Minkowski summand of the zonogon. 

Part (iv) also follows by induction. Using the same argument as for part (ii), note that the only ways to 
have two distinct vertices of the hypercube Tik (of dimension k) project onto the same vertex of the zonogon 
Q is to either have this situation happen for one of the two k — 1 dimensional hypercubes (in which case 
the induction hypothesis would complete the proof), or to have zero translation between the two zonogons, 
which could only happen if ak = bk = 0. 

Part (v) follows by using parts (iii) and (iv) and the definition of a zonogon as the Minkowski sum of 
line segments. In particular, since the difference between two consecutive vertices of the zonogon, Vi^Vi^i, 
for the case k = p, is always given by a single column of the projection matrix (i.e. [aj,bj]'^, for some j), 
then the unique vertices oiTLk that were projecting onto Vi and v^+i, respectively, must be incidence vectors 
that differ in exactly one component, i.e. are adjacent on the hypercube TLk- O 

7.3 Technical Lemmas. 

This section of the Appendix contains a detailed proof for the technical Lemma 4.3 introduced in Section 
4.1.2, which we include below, for convenience. 

Lemma 4.3. When the zonogon has a non-trivial intersection with the band Blu (case [C4]j, the 
convex polygon Ar and the set of points on its right side, r-side(Ar), satisfy the following properties: 

1. r-side(Ar) is the union of two sequences of consecutive vertices (one starting at vq, and one ending at 
Vk ), and possibly an additional vertex, it '■ 

r-side(Ar) = {■mq, ■Wi, . . . , U {vt} U {vr, Vr+i Vk} , for some s < r e {0, . . . ,k}. 

2. With cotan(-, •) given by (25) applied to the (71,72) coordinates, we have that: 

cotan('Us, 'Umin(t,r)) ^ ' whcncvcr t>s 
cotan (■u,nax(t,s), Vr) < whenever t <r. 

Proof of Lemma 4-3. In the following exposition, we will use the same notation as introduced in Section 
4.1.2. Recall that case [C4] on which the lemma is focused corresponds to a nontrivial intersection of the 
zonotope 8 with the horizontal band Blu defined in (30). As suggested in Figure 5 of Section 4.1.2, this 
case can be separated into three subcases, depending on the position of the vertex Vt relative to the band 
Blu, where the index t is defined in (31). Since the proof of all three cases is essentially identical, we will 
focus on the more "complicated" situation, namely when Vt S Blu- The corresponding arguments for the 
other two cases should be straightforward. 

First, recall that Ar is given by (28), i.e. Ar = convdi^o, . . • ,Vk}), where the points Vi are given by 
(29), which results from applying mapping (27) to Vi £ &. From Definition 4.1 of the right side, it can be 
seen that the points of interest to us. namely r-side(Ar ), will be a maximal subset {'^^(i) , '^^(2) : ■ ■ ■ 1 Vi(m) } ^ 
{vo,---,Vk}, satisfying: 

*i(i)=argmax argmin[7 : 7 = (71, 72) £ {i)o, ■ • ■ , ] 

71 72 

= argmax argmax [7 : 7 = (71, 72) e {vq, ■■■,Vk} ] (65) 

71 72 

cotan (■Ui(i), i)j(2)) > cotan {v^2), *»(3)) > • • • > cotan ({>i(„i_i), {'i(„i)) . 
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For the analysis, we will find it useful to define the following two indices: 

s = mm{i G {0, . . . , fc} : O^iv,) > y* ~ U} , f max{i £ {0, . . . , fc} : e2{v,) < y* ~ L} . (66) 

In particular, s is the index of the first vertex of r-side(0) falling inside Blu, a-nd f is the index of the last 
vertex of r-side(0) falling inside Blu- Since we are in the situation when Vt G Blu, it can be seen that 
0<s<t<f<k, and thus, from (31) (the definition of t) and (32) (typical conditions for the right side of 
a zonogon): 

^>...>^>...>^>c>^>...>^>...>^. (67) 

bi bs bt bt+i bf bk 

With this new notation, we proceed to prove the first result in the claim. First, consider all the 
vertices Vi G r-side(8) falling strictly below the band BlUi i-e. satisfying 6*2 [fi] < y* — U. From the 
definition of s, (66), these are exactly dq, . . . , fs-i, and mapping (27) applied to them will yield: Vi = 
( 9i[vi] + c - U, 62[vi] + U ). In other words, any such points will simply be translated by (c • U,U). Simi- 
larly, any points Vi G r-sidc(0) falling strictly above the band Blui i-C- d2[vi] > y* — L, will be translated 
by (c ■ L,L), so that we have: 

v,=v, + {c-U,U), i G {0, . . .,s - 1}, 
^ v., + {c- L,L), ie {f' + l,...,k}, 

which immediately implies, since Vi G r-side(0), that: 

J cotan('Uo, ■"!) > cotan(ui, V2) > ■ ■ ■> cotan('Ms-2, t^j-i) , 
I cotan('Ur+i, Vf+2) > cotan('Ur+2, ■Wf'+s) > ■ • • > cotan('Ufe_i, Vk) ■ 



(69) 



For any vertices inside Blu, i-e- Vi G r-sidc(0) n Blu, mapping (27) will yield: 

Vt = { 0i[vi] - c- 02[vi\ + c- y* , y* ) , i & {s, . . . ,t, . . . ,r}, (70) 

that is, they will be mapped into points with the same 72 coordinates. Furthermore, using (23), it can be 
seen that Vt will have the largest 71 coordinate among all such vt: 

li[it] ~ li[it] = Oi[vt] - ei[vt] - c - {02[vt] - 02 bd) 

(23) / E*=»+i % - c ■ bj^> 0, iis<i<t 

(67) V'-Lj 

E}=t+i aj+c- bj > 0, ift<i< r. 

Furthermore, since the mapping (27) yielding 72 is only a function of 62, and is monotonic non-decreasing 
(strictly monotonic increasing outside the band Blu), vertices t>o, ■ ■ ■ ,Vk G r-sidc(0) will be mapped into 
points {jo, A; G r with non-decreasing 72 coordinates: 

12[vq\ < J2[vi] < - - - < 72['Us-l] < y* = j2[Vs] = • - - = 72 [■&*] = ■ - - = 72 < 72['»>f+l] < ' ' ' < ■y2[Vk]- 

Therefore, combining this fact with (69) and (71), we can conclude that the points Vi satisfying conditions 
(65) are none other than: 

r-sidc(Ar) = {^O,*!, - ■ ■ ,Vs,Vt,Vr,Vr+l,Vk} , 

where the indices s and r are given as: 

dof I max{i G {1, . . . ,s — 1} : cotan (vi-i, vA > cotan({>i, Vt)} 
s = < . 

0, if the above condition is never true, 

(72) 

dof J niin {i G {f + 1, . . . , fc — 1} : cotan (■Wt, ■Si) > cotan ({ji, Vi+i)} 
I A:, if the above condition is never true. 
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This completes the proof of part (i) of the Lemma. We remark that, for the cases when Vt falls strictly 
below Blu or strictly above BlUi one can repeat the exact same reasoning, and immediately argue that the 
same result would hold. 

In order to prove the first claim in part (m), we first recall that, from (72), if s < s — 1, we must have: 

cotan('Us, 'Us+i) < cotan (vs-i-i, Vt) , 
since otherwise, we would have taken s + 1 instead of s in (72). But this immediately implies that: 

cotan(t>s, t's+i) < cotan(t)s+i, Vt) ^ — : — t < — : — r ^ (mediant mcquality) 

72 [Vs+l\ - 72 [Vs\ 72 [Vt\ ~ 7l Fs+lJ 
Jl[Vs+l]-ll[Vs] . ll[Vt]~^l[Vs] (68) a.s+1 , , ,~ ~ s 

which is exactly the first claim in part (m). Thus, the only case to discuss is s = s — 1. Since s > 0, it must 
be that, in this case, there are vertices Vi G r-sidc(0) falling strictly below the band Blu- Therefore, we can 
introduce the following point in 0: 

M = arg max { (01 , 02) e e : 02 = 2/* - [/ } (73) 

Si 

Referring back to Figure 6 in Section 4.2, it can be seen that M represents the point with smallest 02 
coordinate in Blu H r-side(O), and M G [vs~i,Vs]- If we let (0i [Mj, 02[M]) denote the coordinates of M, 
then by applying mapping (27) to A/, the coordinates of the point M G F are: 

M = {9i_[M]+c-U, 02[M] + U) = {ei[M]+c-U, y* ) . (74) 

Furthermore, a similar argument with (71) can be invoked to show that 71 [-A/] < 71 [{>t]. With s = s — 1, we 
then have: 

cotan(t)s, Vt) = — ^ ^— > (smce 72[ftJ = 72[A^J = y > 72[fs-iJ) 

72l^^tJ - 72Fs-lJ 

^ 71 [M] - 7i[vs-i] 



72 [M] - 72[i's-i] 

(68M74) 01 [Af] - Bi{Vs-i\ 
02[Af]-02b,_l] 



= (since M G [t)j_i, v^]) 



which completes the proof of the first claim in part (ii). 

The proof of the second claim in (ii) proceeds in an analogous fashion, by first examining the trivial case 
r > f+1 in (72), and then introducing N argmaxg^ { (0i, 02) G © : 02 = y* — } for the case r = □ 
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